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ABSTRACT 

Extremely red objects, identified in the early Spitzer Space Telescope obser- 
vations of the bright-rimmed globule IC 1396A and photometrically classified as 
Class I protostars Class II T Tauri stars based on their mid-infrared colors, were 
observed spectroscopically at 5.5-38 /im ( Sj)zfeer InfraRed Spectrograph), at the 
22 GHz water maser frequency (NRAO Green Bank Telescope), and in the optical 
(Palomar Hale 5-m), to confirm their nature and further elucidate their proper- 
ties. The sources photometrically identified as Class I, including IC1396A:a, 7, S, 
e, and (, are confirmed as objects dominated by accretion luminosity from dense 
envelopes, with accretion rates 1-10x10~^Mq yr~^ and present stellar masses 
0.1-2 Mq. The Class I sources have extremely red continua, still rising at 38 fim, 
with a deep silicate absorption at 9-11 /im, weaker silicate absorption around 18 
/im, and weak ice features including CO2 at 15.2 /im and H2O at 6 /im. The 
ice/silicate absorption ratio in the envelope is exceptionally low for the IC 1396 A 
protostars, compared to those in nearby star-forming regions, suggesting the en- 
velope chemistry is altered by the radiation field or globule pressure. Only one 
22 GHz water maser was detected in IC 1396A; it is coincident with a faint mid- 
infrared source, offset from near the luminous Class I protostar IC I396A7. The 
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maser source, IC 1396A7fe, has luminosity < O.ILq, the first H2O maser from 
such a low-luminosity object. Two near-infrared H2 knots on opposite sides of 
IC 1396A:7 reveal a jet, with axis clearly distinct from the H2O maser of IC 
1396A:76. 

The objects photometrically classified as Class II, including IC1396A:/9, 6, 
2MASSJ 21364964+5722270, 2MASSJ 21362507-^5727502, LkHa 349c, Tr 37 
11-2146 and Tr 37 11-2037, are confirmed as stars with warm, luminous disks, 
with a silicate emission feature at 9-11 /im, and bright Ha emission, so they 
are young, disk-bearing, classical T Tauri stars. The disk properties change 
significantly with source luminosity: low-mass (G-K) stars have prominent 9-11 
emission features due to amorphous silicates while higher-mass (A-F) stars have 
weaker features requiring abundant crystalline silicates. A mineralogical model 
that fits the wide and low-amplitude silicate feature of IC1396A:6' requires small 
grains of crystalline olivine (11.3 /xm peak) and another material to to explain 
the its 9.1 /zm peak; reasonable fits are obtained with a phyllosilicate, quartz, or 
relatively large (> 10 /xm) amorphous olivine grains. 

The distribution of Class I sources is concentrated within the molecular glob- 
ule, while the Class II sources are more widely scattered. Combined with the 
spectral results, this suggests two phases of star formation, the first (4 Myr ago) 
leading to the widespread Class II sources and the central O star of IC 1396, and 
the second (< 1 Myr ago) occurring within the globule. The recent phase was 
likely triggered by the wind and radiation of the central O star of the IC 1396 
H II region. 

Subject headings: young stellar objects; ISM: individual (IC 1396A); ISM: glob- 
ules; infrared: stars; planetary systems: protoplanetary disks 



Introduction 



Spitzer has opened widely the mid-infrared (3.6-24 /im) window for studies of star 
formation. In the earliest observations with the observatory, as in many subsequent ones, 
previously unknown, mid-infrare d-bright sources h ave been discovered in star-forming re- 
gions. In this paper we follow up (iReach et al.ll2004l . Paper I) on the mid-infrared sources in 
the bright-rimmed globule IC 1396A, known as the Elephant Trunk Nebula, a dense globule 
in the large, nearby (750 pc) H II region IC 1396, which is is excited by the 4 Myr-old 06 
star HD 206267 (see IWeikard et alJll996r). Based on cornparis on to measured colors to those 
of young stellar objects in Taurus (IKenyon fc Hartmannlll995[ ). we suggested 10 sources may 
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be the very early Class I stage, which is thought to only last 10^ yr (jStahler &: Pallal 12004 ). 
indicating that the globule is a site of very recent star formation. 

Detecting new protostars in a wide variety of star forming regions is an important 
expansion of previous work that by necessity has concentrated on the very closest star forming 
regions such as Taurus and Orion. While these nearby regions provide excellent opportunities 
to observe the widest range of star formation (from brown dwarfs to massive stars), they 
are necessarily parochial. Contrary to a long-standing common belief that the Sun formed 
in a quiescent region like Taurus, meteoritic evidence indicates the Sun likely formed near 
a high-mass star, whose H II region, winds, a nd supernova explosion influence the evol ution 
and composition of the protoplanetary disk (IHester et al.ll2004j : iTachibana et al.ll2006l ). To 
search for present-day star-forming regions similar to those in which the Sun formed, we 
must consider massive star forming regions. I n such regions, stars can b e formed by radiative 
driven implosion of moderately dense cores Jhefloch fc Lazarefjll994h. and planet-building 
disks are exposed to strong ultraviolet radiation (jo'Dell fc Wen 1994 ). A major goal of 
future star formation studies should be to determine the effects of massive stars both in 
triggering star formation and shaping disks around young stellar objects. 

In this paper, we present follow-up observations to examine the nature of the sources in 
IC 1396A in more detail. New Spitzer mid-infrared spectra demonstrate that the photomet- 
rically extremely red objects are indeed Class I protostars, with dense, highly-optically-thick 
envelopes shrouding a moderately luminous core inside a dust photosphere. We show that 
at least one of the extremely red objects, IC 1396 A:7, powers a molecular outf low with H2O 
mase r emission, characteristic of a very young (~ 10"^ yr) Class protostars (IFuruya et al. 
2OOII ). The photometrically-classified Class II objects are compared to classical T Tauri 
stars, based on Ha emission and mid-infrared disk spectra. 



2. Observations 
2.1. Spitzer Space Telescope 



Mid-infrared spectra were obtained using the InfraRed Spectrograph (IRS; iHouck et al. 



20041 ) in 2004 and 2005 November. Table [T] lists the sources and IRS observing modes. The 
2004 observations was adversely affected by poor pointing; the 2005 observations repeated 
those of the first set that had offsets between actual and expected pointing larger than half 
the slit width in the cross-slit direction. Spectra were taken with both the short- and long- 
wavelength, low-resolution mode for all sources. The brightest sources were also observed 
with the high-resolution mode. Observations were performed in groups ('clusters'), and a 
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single bright, isolated mid-infrared source (whose coordinates were measured from the IRAC 
images taken in 2003 December) was used as a peak-up source for each group. Mid-infrared 
source brightnesses and positions were verified using the IRAC and MIPS images from Paper 
I. The magnitudes are reported in Table [2J With respect to the values in Paper I, significant 
changes include the following. Source IC 1396A:A is significantly fainter in the IRAC 3.6 
and 4.5 /xm channels, because in Paper I it was confused with another nearby source that is 
brighter in those channels (but far fainter at 24 /im). The changed colors move the source 
in the color-color diagram (Fig. 3 of Paper I) from an unlabeled region into the Class 
1 region. The position of source IC 1396A:k was corrected from the typographical error 
(copied declination of A). 

For the basic reductions, the on-slit nods wer e subtracted from each other, and spectra 



were extracted using SMART (IHigdon et al.ll2004l ). In many cases, the nod-subtracted spec- 



tra were contaminated by other sources or bright structures in the nebula; for these, a local 
background subtraction was performed on each nodded image. The extracted spectra from 
each nod position were averaged, and their uncertainties combined using the root-sum-square 
of each extraction's uncertainty. 

For IC 1396A:a;, 7, 6, and e, the low-resolution [R = A/AA ~90) 2-dimensional Ba- 
sic Calibrated Data (BCD) spectral images produced by the SSC pipeline versions S15.3.0 
and S17.2.0 were used as the starting point of a customized reduction. First, the spectral 
and spatial dimensions were orthogonalized, and then the 2-dimensional images of the two 
nodding positions were subtracted in order to remove extended background emission. In 
some cases the background emission is highly structured or multiple sources are present 
in the slit, and nodding pairs could not be subtracted. Here, the background observed in 
nearby telescope pointings was subtracted and any residual emission was removed using the 
background determined from neighboring columns in the source extraction process. Some 
background emission lines (especially at 6.2, 17.0, and 18.7 fim) were not removed completely 
or were over-corrected for, which was taken into account in the interpretation of the YSO 
spectra. The background subtraction process also removes most of the effects of pixels with 
deviating dark currents ("hot pixels"). Any remaining bad pixels were replaced by values 
from a two-dimensional interpolation. Subsequently a fixed width (4 pixels) extraction was 
performed and the 1-dimensional spectra were then averaged. Then the spectra were divided 
by spectra of the standard star HR 2194 reduced in the same way in order to correct for 
wavelength-dependent slit losses. Spectral features in HR 2194 (AO V) were divided out 
using the photospheric model of Decin et al. (2004). Some residual standing wave structure 
with a period of ~1 /im observed above 20 fim was not corrected for. Finally, the Spitzer/IRS 
modules and ground-based spectra were multiplied along the flux scale in order to match 
Spitzer/IRAC 3.6 and 8.0 /im fluxes, using the appropriate filter profiles, and to provide 
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smooth transitions between the IRS modules. 



2.2. Keck/NIRSPEC L-band Spectra 

The Spitzer spectra of IC 1396A a and 7 are complemented by ground-based L-band 
spectra obtained wit h the NIRSPEC spectrometer at the Keck II telescope at Mauna Kea 



(iMcLean et al.lll998l ) at a resolving power of R = A/AA=2000. This provides an independent 
measure of the solid H2O column density through the strong 0-H stretching mode at 3.0 fim. 
The data were obtained on UT 24 July 2004 in two separate grating settings covering the 
2.70-3.54 and 3.34-4.19 /xm spectral ranges. The telescope was nodded on the sky in an ABBA 
pattern for 12 and 8 minutes for the short and long-wavelength settings respectively. The 
spectra were reduced in a way standard for ground-based long-slit spectra, using the nearby 
main sequence star HR 8208 (FOV) for telluric line correction. The final signal-to-noise values 
are 30 and 50 for the short and long wavelength sections, respectively. Finally, the ground- 
based spectra were multiplied along the flux scale in order to match the Spitzer/IRAC 3.6 
/xm photometry, using the appropriate filter profile. 



2.3. Palomar Hale Telescope: Optical spectra 



Visible spectra were obtained using the Double Spectrograph (jOke fc Gunru Il982l ) in 
August 2006. Visible counterparts were identified for the Class II sources using finder charts 
from the Palomar Observatory Sky Survey (POSS), 2nd-generation Digitized Sky Survey; 
these sources were identified easily and placed in the long slit with 1" width. For objects 
not visible on the POSS, offsets were performed from a nearby star, using the coordinates 
from the IRAC images relative to the star in the POSS. All spectra were obtained in the 
low-resolution mode, with a resolving power 2100 (2500) in the blue (red) portion of the 
spectra. A dichroic with transition at 5500 A was used to separate light to the blue and 
red gratings and cameras. The blue spectrum was centered at 4400 A, with useful range 
2700-5500 A. The red spectrum was centered at 6800 A, with useful range 5500-7900 A. The 
spectra were analyzed using IRAF and then normalized in IDL. For objects not visible on 
the POSS, their spectra were extracted using a reference position on the slit from a bright 
sources. 
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2.4. Green Bank Telescope 

Spectroscopic observations to search for H2O masers from fourteen mid- infrared sources 
in IC 1396A were performed using the National Radio Astronomy Observatory (NRAO)QlOO- 
m Robert C. Byrd Green Bank Telescope. The targets were observed in multiple sessions, 
to allow for the possible strong variability of protostellar masers on month timescales. The 
observing sessions were on 2006 Aug 7 and 2006 Sep 27. 

The 18-22 GHz receiver was tuned to 22.23508 GHz, and the spectroscopic observations 
were taken in two methods. First, each source was observed by nodding between the two 
beams of the receiver; this efficient technique always has the source within one of the beams, 
but it only works the 'off' beam located 3/ in azimuth is free of emission. Second, each source 
was observed by position-switching to a blank (in radio and mid-infrared emission) field. 
Both the nods and position-switched spectra were analyzed using GBTIDL. The well-known 
maser source S Per was observed to validate the observing technique, and it was cleanly 
detected with an antenna temperature 140 K, consistent wit h the previously-measured gain 



of 1.8 K/ Jy and the previously-measured S Per fiux of 76 Jy (jVlemmings fc van Langevelde 



20051). 



In all cases, the nodded and position-switched spectra agreed, and for all spectra except 
those near IC 1396A:7 the spectra are consistent with noise. Upper limits for undetected 
sources were derived as follows. Each 1-min-per-beam nodded spectrum yielded an rms 
noise of 0.24 K antenna temperature with 0.01 km s~^ autocorrelator bins after banning 
smoothing. Rebinning to 0.1 (1) km s^^ bins, the rms antenna temperature was 0.11 (0.04) 
K. Therefore, we place a 3-cr upper limit on the H2O masers for all observed sources (except 
detections noted below) of 0.18 Jy peak (in 0.1 km bins) and 0.07 Jy km s~^ fiux (in 1 
km s~^ bins). For sources IC 1396A:7, a, and 6, observations were 8 times longer yielding 
an rms of 0.028 K, and a 3-0" upper limits of 0.05 Jy peak and 0.02 Jy km s~^ fiux. 

Figure [T] shows the H^O maser spectra of IC1396A:7. This source had been previously 



detected by IValdettaro et al.l (j2005[ ) using the Medicina telescope with 1.9' beam. Relative 
to the peak fiux of 33.7 Jy km s~^ in 2004 Dec as observed at Medicina, the maser fiux 
decreased by more than a factor of 10^ in 2006 Aug, the first GBT epoch. By the last GBT 
epoch in 2007 Jan, the maser brightness had increased, but still only to a factor of 120 less 
than the fiux in 2004 Dec. Continued monitoring at Medicina at 7 different epochs during 
2007 and in 2008 Apr yielded upper limits of 1-2 Jy, confirming that the maser has remained 



^The National Radio Astronomy Observatory is a facility of the National Science Foundation operated 
under cooperative agreement by Associated Universities, Inc. 
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faint since 2006. 



The GBT observations with 33" beam show the maser is very close to IC1396 A:7. 
A cross pattern of observations through the source shows the offset of the centroid of 22 
GHz maser emi ssion from the in fr ared source is 11" ± 11". High-angular- resolution VLA 
observations by IValdettaro et al.l (120081 ) showed that the H2O maser is actually 7" N of 
IC1396A:7. The maser spot is coincident with another, fainter mid-infrared source. The 
VLA observations were taken in 2005 Feb, when the maser was very bright (73 Jy), and the 
beam size was 0.4" so the shift of 7" is highly significant. 
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Fig. 1. — GBT 22 GHz H2O mascr spectra of IC1396A:7 on 4 dates (as labeled). The spectra 
are offset vertically by 0.2 Jy on each successive date, for clarity. 
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2.5. 2MASS 



The 2 Micron All-Sky Survey (2MASS; Iskrutskie et~al1l2006h Extended Mission @ pro- 



vides additional products to the 2MASS survey. For a few regions of the sky, including IC 
1396A, 6 time exposures have been obtained. As a longer exposure time extends the com- 
pleteness level by about 1 mag, there is a drawback regarding the photometric saturation 
which occurs for fainter sources. Consequently, photometry from the main 2MASS survey 
must be preferred for all saturated sources in the 2MASS 6X catalog (i.e. when the flag 
rd_f Ig equals 1). 



We built the extinction map of IC 1396A following the method described in lCambresy et al 



(120021) using the 2MASS H - color excess of 32000 stars in a 40' x 40' field around the 
nebula. Sources detected at 24 /im have been removed from the catalog since they are intrin- 
sically red objects for which the color is not related to the cloud extinction. The extinction 
is obtained following the expression: 

where {Ah/ Ay - AkJAy)'^ = 15.87 (Rieke & Lebofsky 1985) and Z^^i = -2.5 mag. 
The zero point for the extinction calibration, Zco\, corresponds to an intrinsic star color of 
H — Ks = 0.16 mag which is measured outside the nebula. Any diffuse extinction in the 
field is therefore removed by construction. The extinction is obtained using the median of 
the 10 stars included in each cell. It produces a map of IC 1396A with a resolution of about 
1 arcmin, with an average visual extinction of about 5 mag. The extinction toward some 
specific Class II objects, which will be used later in this paper are Ay = 5.4, 3.8, 3.5, 4.5, and 
3.5, for IC 1396A:^, P, Tr37 11-2037, Tr37 11-2146, and 21364964+5722270, respectively. 

The extinction map clearly shows the central hole associated with LkHa 349 in the 
higher density round head of the globule. The visual extinction in the hole is 2.8 mag 
whereas it is ~ 6 mag in the north side of the head and ~ 10 mag for the south side 
with a maximum of 16.3 mag at the position 21h36m59s +57d30m00s (J2000). Beside, the 
correlation of the dust extinction and the emission in the Spitzer images exhibits a region 
with a deficit of infrared emission that corresponds to the ~ 5' size structure at 21h36m 
+57d23m, LDN 1099. This deficit can be explained assuming LDN 1099 is disconnected 
from IC 1396A and located at a larger distance of the 06 heating star. This is confirmed by 
the examination of the optical images which suggests IC 1396A is illuminated from behind 
whereas L1099 would be illuminated on its front size. Assuming the system composed by 



^ ht t p : / / www . ipac .caltech . edu / 2niass /releases/ allsky /doc/ explsup . html 
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the 06 star, at 750 pc, and IC 1396A is seen with an angle of 45 degrees (0 degree meaning 
the star and the globule are at the same distance), we find the main globule IC 1396A is 4 
pc closer to us than the ionizing star and that L1099 is 9 pc farther. This is confirmed by 
the difference of radial velocities fo r the two regions wh ich are -8 km s~^ and km s~^ for 
IC 1396 A and L1099, respectively JWeikard et aljll996h . 



2.6. Palomar Hale Telescope: Near-infrared images 



IC 1396A was imaged using the Wide-Field Infrared Camera (WIRC: IWilson et al.ll2003l ) 
on the 200" Hale telescope on 2004 Aug 25-26. The 8.7' field-of-view chip was steered toward 
a series of 22 positions alternating between a gradual scan across the globule and a nearby 
empty reference field. The reference field images were combined to generate a sky fiat, which 
was then removed from each image. While observing, a small dither offset was added to 
each position to ensure that stars in the 11 reference field observations were removed during 
creation of the sky fiat by robust averaging. The astrometry for each image was measured 
by comparing to the 2MASS point source catalog. Then images were combined into mosaics 
for each filter. 

The absolute calibration of the mosaics was measured by correlating photometry in 
image units to the corresponding entries in the 2MASS catalog. The cross-calibration is 
accurate to 5%. The location of each of the sources discussed in this paper was inspected 
in each of the near-infrared images. Some sources that were marginally detected by 2MASS 
(e.g. a, 7, 7], in J band) or suffered from blending with nearby sources {(3) were clearly 
detected in the Palomar images (0.9" FWHM seeing). Table [2] (columns J and K) lists the 
magnitudes determined from the Palomar observations. 

Figure [2] shows the near-infrared image. All of the bright Class I sources for which we 
obtained mid-infrared spectra are detected in this image. Two exceptionally red sources, a 
and 7, may be marginally resolved in the Palomar J-band image, with sizes ~ 1.5" (1100 
AU). In all the Palomar near- infrared images, the protostars are dominated by continuum, 
since their fluxes in line (Pa/5 and H2 2.12 /xm) and continuum (narrow 2.17 fim and broad 
J) filters are similar. The resolved J-band emission is likely to be scattered light. 



Fig. 2. — Near-infrared image of IC 1396A from Palomar, combining the 2.2 /im narrow con- 
tinuum (red), H2 (green), and J-band (blue). Since the H2 and 2.2 /xm continuum wavelengths 
are similar, point sources with spectral energy distribution rising to longer wavelengths ap- 
pear red-|-green=yellow. The Class I protostars for which Spitzer spectra are presented in 
this paper are labeled, as are the LkHa 349 stars and the Class II star Tr 37 11-2146. The 
mid-infrared sources are generally 'yellow' in this representation, but they do not stand out 
based on near-infrared colors alone. The vast majority of near-infrared 'yellow' sources, 
which preferentially fall within the globule, are faint or not detected at 8 /im with IRAC 
(too faint for mid- infrared spectroscopy with IRS); they are probably extincted background 
stars. 
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3. Class I Sources 

3.1. Properties of the mid- infrared spectra 

Figure [3] shows the spectra of the bright Class I protostars a, 7, 6, e, i], and A. Table O 
lists several properties derived from the mid-infrared spectra including colors and silicate 
feature amplitude. For the Class I sources, the silicate features at 10 and 18 yum are in 
absorption hence AF/F is negative; conversely for most Class II sources the silicate features 
are in emission. 
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Fig. 3. — Spitzer/IRS spectra of Class I sources IC1396A:A, r], a, 7, 5, and e. Each spectrum 
is in uLi, units, divided by its value at 8 /im and offset vertically by a factor of 3.7 for 
clarity. Spitzer and Palomar photometry are over-plotted as asterisks. For sources a and 7, 
IRS high- resolution spectra and Keck L-band spectra are included. For source a, the IRS 
high-resolution spectrum continuum does not agree with that in the low resolution-spectrum 
from 15-20 /xm; the low-resolution continuum is believed to be more accurate (and agrees 
with the MIPS 24 /xm photometry) . Narrow spikes in the spectra are due to rogue pixels and 
should be ignored. The positive feature in A and a and negative dip in 77 at 11.2 /im are due 
to incomplete subtraction of a bright, nebular PAH feature; the wavelengths of nebular hues 
that are incompletely removed in some spectra are indicated by dotted lines at 6.2 (PAH), 
11.2 (PAH), 12.8 ([Nell]), and 18.7 ([S III]) pim. 
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3.2. Ice absorption features 



The Spitzer spectra of IC 1396A a, 7, 6, and e show the well-known ice absorption 
features between 5 and 8 fim, and at 15.2 fim (Fig. [3]). The latter is due to solid CO2, 
while the 5-8 fim ab sorption is caused by the overlapping bands of many solid state species 
(iBoogert et al.ll2008l ). At the long- wavelength side, the 5-8 /im absorption overlaps with 
the edge of the strong 9.7 fim band of silicates, which complicates the determination of an 
overall continuum level for the study of the ice absorption features. A best-effort global 
continuum was determined, using all sp ectra and pho t ometr y available between 1 and 40 
/xm, similar to the method described in iBoogert et al.l (120081 ). It is kept in mind that the 
optical depth of the broad absorption overlying the entire 5-8 /im range shown in Fig. H] is 
relatively uncertain, but that of individual, narrower, features within it is not. 

The bending mode of solid H2O peaks at 6.0 /im, but is a prominent contributor to the 
entire 5-8 fim absorption region. In order to determine its contribution, the H2O column 
density is obtained independently from the 3.0 /xm H2O stretching mode, when available 
(sources a and 7) or from the 13 fim libration band (Table [6]). Clearly a significant fraction 
of the 5-8 fim absorption is not due to the bending mode of solid II2O (Fig. H]). The excess 
absorption between 5.7 and 6.6 fim is likely real, and is caused by the overlapping bands of 
solid H2CO (5.7 /im), HCOOH (5.8 /im), NH3 (6.12 /xm), and possibly ions (e.g. HCOQ-) 
and PAH species near 6.3 fim. 

At longer wavelengths, a distinct band centered at 6.9 /im is detected toward all sources. 
This "6.85 /im band" is commonly observed toward low- and high-mass YSOs, as well as 
background stars. Its profile is known to vary considerably and has been decomposed in short 



and long-wavelengths components (called C3 and C4 in iKeane et al.l I2OOII : IBoogert et al. 



20081 ) centered on 6.755 and 6.943 /im, respectively. The 6.85 /im bands toward the IC 
1396A protostars are shifted remarkably f ar to longer wavelengths, i.e. C4 is much stronger 
than C3. Previously (IBoogert et al.ll2008l ) it was found that such sources tend to have low 
solid H2O abundances, expressed in the H2O column density scaled to the peak optical depth 
of the observed 9.7 /zm band of sihcates. 

Indeed, the IC 1396A protostars have very low A^(H20)/r9.7 ratios as well (Fig. [5]). The 
origin of the components of the 6.85 /im band is not fully established, but the correlation 
with dust temperature tracers, such as shown in Fig. 0, suggests its carrier(s) are affected by 
ice processing. The NH^ ion has been considered a reasonable candida te, as upon heating 



its b and shifts to longer wavelengths and its band strength increases (ISchutte fc Khanna 
2OO3I ). It is also more refractory than H2O ice. Regardless of its origin, it appears that the 
ices toward the IC1396A protostars are thermally, and possibly energetically processed. All 
4 studied YSOs show these unusual ice characteristics, which is distinct from most sources in 
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Fig. 4. — Optical depth spectra (dots) of the Class I objects in IC 1396A showing 5-8 /xm 
absorption bands. The previously studied spectrum of the low mass YSQ HH 4 6 IRS, repre- 
sentative of the large sample of low mass YSOs studied in lBoogert et al.l (120081 ) . is included 
at the bottom. The spectra are normalized to the peak optical depth by multiplication with 
the factor indicated on the left, and then shifted along the vertical axis for clarity. For each 
source the contribution of the bending mod e of pure solid H2O is shown by a laboratory 
spectrum (smooth line; iHudgins et al.lll993l ). Note that in all IC 1396A YSOs the distinct 
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the la rge sample of low mass YSOs in nearby star-forming regions studied by lBoogert et al. 
( 120081 ) ■ This suggests that the ices toward IC 1396A are processed by radiation external to 
the YSOs, likely the harsh radiation from the nearby O star. 

Finally the solid CQ2 abundances relative to H 2O are between 14-34% (Table [6]), which 
is typical for most YSOs (jPontoppidan et al.ll2008l ). Solid CH3OH was looked for at 9.7 and 
3.53 fxm as well, but is less than 5% with respect to H2O. 



3.3. Models for the Class I protostars 



We now attempt to derive the physical properties of the Class 1 protostars from the 
observed spectra and photometry. To this end, we first create a simple 'toy' model of a 
Class I object and use it to derive some fiducial parameters for the envelope. The simple 
model is a combination of modified blackbodies, wi th outer layers a bsorbing inner ones such 
as was applied to protostars in the Trifid nebula (IRho et al.ll2006l ). Then we compare the 



obser vations to more sophisticated radiative transfer models developed by (jWhitney et al. 
20031 ). which allow us to estimate the properties of the central stars and their accretion 



rate. While the theoretical models are more physically detailed, they don't fit the spectra 
as accurately as the simple models, which we present here as an empirical description of the 
data. 



3.3.1. Simple model 

The spectral energy distributions of the Class 1 sources appear to be the sum of two 
primary color temperatures, with the warmer component having Ti ^ 420 K and a colder 
component having T2 ^ 120 K. To explain the 70 /im emission, a colder component is 
required, with T3 ~ 40K. As a first, empirical, description of the spectra, we will model each 
protostar as the sum of three nested components, with the higher temperature ones residing 
closer to the center and being absorbed by the colder ones: 

F°^' = Fje-("^+"») + Fje-"3 + Fj^. (2) 

This simple model neglects direct emission from a potential central star, as well as neglecting 
details of radiative transfer. Neglecting the central star is entirely justified based on the lack 
of optical counterpart and the great optical depth toward the center. Neglecting radiative 
transfer is not justified, so one must bear in mind that the empirical models here are only 
descriptions of the data and not necessarily self-consistent. (More sophisticated models are 
described in the next section.) The optical depth of the outer layers was normalized to 
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Fig. 5. — Correlation of the ratio of integrated optical depths of the components C4 and 
C3 of the 6.85 /im bands with a r aeasure of the solid H^O abundance for YSOs studied 
in the c2d ic e survey (filled cir cles; iBoogert et al.ll2008l ). ISO-observed massive protostars 
(open circles; iKeane et al.ll200ll ) and our results for IC 1396A (triangles). All 4 YSOs in IC 
1396A have strongly red-shifted 6.85 fim bands, i.e. large C4/C3 ratios and very low solid 
H2O abundances. This likely implies a high degree of thermal processing of the ices in the 
IC 1396 A core by the nearby 0-star. The sources with properties overlapping with the IC 
1396A protostars (C4/C3/,1.2 and H2O/r9,7j0.25) from previous studies are (from high to 
low C4/C3): Mon R2 IRS3, IRAS 03301+3057, S140 IRSl, DC Tau B, and W3 IRS5. 
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match the observed depth of the 9-11 /im sihcate feature, and the wavelengt h-dependence 



of the opacity was taken from the interstellar dust model of iDraine fc Lil (1200 ll ). The coldest 
component (3) was fixed at 40 K and was normalized to the difference between the observed 
70 yum flux and the model for components 1 and 2 (the inner and outer envelope, respectively). 

Figure [6] shows the observed spectra and empirical model fits. The color temperatures 
and ratio of the components 1 and 2 were adjusted to match the spectrum from 5-40 fim. 
We found it possible to match all the spectra with similar values of Ti and T2, so we fixed 
them at 420 and 120 K, respectively. Specifying the 22 /im flux, F22, the ratio of fluxes at 
22 to 8 fim, F22/FS, and the silicate feature optical depth, Tsu was sufficient to constrain 
components 1 and 2. Greybody dust grains in region (2) heated by radiation from region 
(1) with luminosity Li will be at temperature 



Combining the outer envelope (region 2) optical depth and size, and assuming dust/gas 
abundance as in the ISM, we then derive the density, n2, in region 2, which is of order 10^ 
cm~^. If the protostar is in steady accretion, then its density should follow 

^2 = . (4) 



[e.g. 



Terebey. Shu, fc CassenI Il984j ) from which we derive mass accretion rates of order 
lO^^M^^ Mq yr~^. The accretion luminosity of the outer envelope alone (i.e. for ini- 
tially distant material that only moves R2) is only of order 0.1 M^,Lq, less than the observed 
luminosity of the outer core {L2 ~ ^Lq). This shows that material must be accreting further 
in than R2. If the accretion luminosity arises at a radius Race, and there is a steady flow of 
matter inward, then 

J-'core D ' ^ ' 

^acc 

which yields accretion radii R 

acc 

30M, AU. 

The coldest component (3), or outer envelope, was fixed at 40 K and was normalized 
to the difference between the observed 70 /im flux and the model for components 1 and 2. 
including for relevant component i the visible extinction, Ay, the distance from the center, 
Ri, the volume density, n^, and the luminosity, Lj. The opacity of the outer envelope was 
small, so the silicate absorption can be identified with component 2, the envelope. The 
volume densit y of the outer envelope is onl y an order of magnitude higher than that of 



a dense core (IMyers. Linke. &: BensonI Il983l ). so it represents the outermost layer of the 



protostar that merges with the surrounding globule gas. 
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Table [7] lists the results of the model fits, In this empirical model, a collapsed object 
of radius < 30/?© has formed, and it resides inside a spherical, optically thick region of size 
~ 2 AU. The mass of the star cannot be readily constrained without further guidance from 
theoretical models. We assume masses of order 1 Mq for the calculations in this section. 
Changing to larger or smaller masses has no effect on the luminosities (which are tied to the 
observed fluxes), but mostly changes the size of the central object (which is hidden within 
the optically thick core). 

The luminosity of sources a and 5 arises primarily from components (1) and (2), which 
is consistent with these objects harboring a central protostar in the Class I phase. The 
luminosity of component (3), the outer envelope, is highest for 7 and e. The key measurement 
is the 70 /xm flux: the ratio of 70/24 /im flux is 22, 2.6, 12, 36 for 7, a, 5, and e, respectively. 
The 70 yum fluxes are uncertain due to the bright nebular emission and large beam; however 
they are sufficiently accurate to reveal differences between the sources. (Sources b and e are 
not cleanly resolved from each other at 70 /im; the peak pixel brightnesses are comparable 
[with e being ~ 7% brighter], so we attributed half the 70 /im flux to each source.) The 
origin of the 70 /im flux from 7 and e may be different. For 7, the presence of an H2O maser 
spot 7" from the infrared source (see §5.2.21) indicates a Class core may contribute some 
of the 70 /xm flux we attributed to 7. For e, the high 70/24 /xm flux (and correspondingly, 
the high luminosity of component (3) in the simple model) is due to the source being more 
envelope-dominated, since the silicate feature is extremely deep for this source. 



We can now compare this empirical, two-component description of the spectral 



e nergy 



distribution to a simple theoretical model for the protostar. Following iStahler fc Fallal (12004J . 
hereafter SP), we expect a dust photosphere, located where the protostar is optically thick 
to its emergent radiation. This 'photosphere' will span a range of radii due to the rapid 
wavelength dependence of the dust opacity. To determine whether the mid-infrared emis- 
sion can be associated with such a photosphere, we solve SP equations 11.12a,b using the 
luminosity 6 L© (like source a) and an assumed mass 1 Mq. The photospheric tempera- 
ture is in the range 500 - 120 K for a mass accretion rate in the range 0.15 - 4xlO~^M0, 
and a photospheric radius in the range 1.4-120 AU, respectively. These temperatures span 
the 'core' and 'envelope' color temperatures required to explain the 5-40 /xm emission. We 
furthermore note that the 'envelope', to which we ascribed the 9-11 /im silicate absorption, 
has, by definition, an optical depth greater than 1 (but not greater than 10) at 10 /xm, 
so by the simplest definitions it must be considered part of the photosphere. That we re- 
quired two color temperatures to match the spectral energy distribution, rather than a single 
blackbody, reflects only that the wavelength-dependence of the opacity leads to an extended 
photosphere. By ascribing the 9-11 [xm silicate absorption to the cooler component, we have 
implicitly assumed a negative temperature gradient, as anticipated for a source powered by 
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Fig. 6. — Empirical model fits to the Class I protostar spectra. For each source a sum (thick 
solid line) of three blackbodies (shown individually as dashed lines), with the warmest one 
absorbed by the cooler ones, matches the observed IRS spectra (thin solid line) and MIPS 
photometry (asterisk). 
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Fig. 7. — Cartoon illustrating the structure of the Class I protostars. The mid-infrared 
emission detected by Spitzer is attributed to the 'dust photosphere' spanning 1-100 AU and 
surrounded by a cooler envelope. 
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accretion luminosity rather than being heated by starhght and cosmic rays from the outside. 
Figure [7] summarizes the protostar parameters derived in this section. The agreement of 
the accretion rate for a simple photosphere approximation with that derived in the previous 
paragraph from the envelope density, and the plausibility of the protostellar radius required 
for accretion to provide the observed luminosity, suggest the simple protostar model pre- 
sented here is at least roughly consistent with growing star surrounded by a thick envelope 
and powered by accretion. 



3.3.2. Association with Whitney protostar models 



Detail ed radiative transfer m odels of protostars with realistic geometries have been de- 
vel oped bv IWhitnev et al.l (120031 ) and calculated over a wide range of physical parameters 
by iRobitaille et al.l (120061 ). To find the Whitney model that best matches the Class I source 
spectra, we first tabulated, for the observed spectra for sources a, 7, 5, e, C, and 77 at a set 
of 14 wavelengths (1.22, 1.63, 2.19, 3.56, 4.51, 5.58, 7.65, 9.95, 12.93, 17.72, 24.28, 29.95, 
35.06, 71.42 /im). These wavelengths were chosen to sample the SED and largely avoid ice 
features; however, the 12.93 yum feat ure possibly contain s a strong H2O libration feature 
whic h is not included in th e models (iBoogert et al.ll2008l ). For each of the 2 x 10^ models 
from Robitaille et al. (^20061 ). we used the model fluxes integrated over the J, H, K, IRAC 3.6, 
IRAC 4.5, and MIPS ToiTfllterl, and individual model fiuxe. at the remaining Spitzer IRS 
wavelengths provided separately by T. Robitaille (personal communication). A goodness- 
of-fit was calculated from the difference between the observed and model tabulated fluxes, 
including for each protostar model 10 possible inclination angles (87.1° [edge-on], 81.4, 75.5, 
69.5 ,63.3, 56.6, 49.5, 41.4, 31.8, 18.2 [pole-on]), variable (1 < Ay < 8 mag) line-of-sight 
redde ning applied with the Milky Way Ry = 3.1 extinction curve (jWeingartner fc Draine 
2OOII ). and a normalization factor within the range ±20% relative to that appropriate for a 
distance of 750 pc. This procedure is very similar to the one available from Robitaille's pro- 
tostar model fitter^, and we have verified for source IC1396Aa that the model identifications 
are in agreement. The goodness-of-fit (x^) was often poor, potentially due to coarseness 
and inapplicability of the model grid, and possible source multiplicity (IC1396A:7 discussed 
below). In many instances, particularly for sources IC 1396A: 6 and e, the large is caused 
by a steeply rising MIPS 70 /im flux that was not reproduced by the models, although the 
rest of the SED was fitted very well. These sources are located in the southwestern rim of IC 
1396A and have large fitted extinctions of Ay ~ 7-8 mag. We suspect that outer envelope 
heating by the external UV radiation field that is not accounted for in the Robitaille models 



•^available at 



http: / / caravan, astro . wise .edu / protostars / repository / index. php 
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is the cause for the enhanced MIPS 70m fluxes. 

The 10 models with the best were inspected in detail for each source. Table [S] shows 
the most important parameters for the best-fitting model for each source. Figure [S] compares 
the radiative transfer model to the observed spectrum of IC1396A:q;. At wavelengths longer 
than 10 fim , the Class I protostar models are dominated by the 'envelope,' while at 5-8 
fim the flux is mostly from the 'disk' in these models. The 9-11 fim silicate feature is due 
both to absorption of the 'disk' by the envelope and radiative transfer through the envelope 
(whose temperature decreases radially outward). Scattered light contributes significantly 
from ~ 2-4 /xm, and emergent flux from the central star becomes dominant at wavelengths 
less than 2 /xm. It is immediately evident that these radiative transfer models are more 
realistic and contain more of the appropriate physics than the simple model in the previous 
section: no direct flux from the star would be detectable in the simple, spherical model; nor 
is the radiative transfer through the envelope taken into account. The detection of these 
sources in the near-inf rared demonstrates that the envelope and/or disk must be flattened 



(or inhomegenous; see llndebetouw et al.l 120061 ). in order to allow scattered light from the 
star to emerge. 

The properties of the central stars are not well constrained, but the model fits clearly 
prefer the youngest available stellar photosphere models, with ages or order 10^ yr for a, 
6 and e, 10^ yr for 7 and 1], and 10^ yr for (. The evident progression of SED shapes in 
Figure [3] agrees with the model fitting results, suggesting that the SED shape is sensitive to 
stellar age. The central star only contributes significantly to the K-band and 3.6 nm IRAC 
band, and it is highly extincted. Taken at face value, the masses are similar, at 0.2-0.3 Mq. 
Since the accretion rates are still high, these may not be the final masses of the young stars 
that will evolve toward the main sequence. If the sources remain at this accretion rate for 
a nominal Class I lifetime of ~ 10^ yr, the added mass would be in the range 0.2-0.9 Mq, 
suggesting the final masses of the stars would be in the range 0.6-1.1 Mq. 



3. 3. 3. Comparing the empirical and radiative transfer model results 

For sources IC1396A:a and e, the total luminosities obtained from the empirical (Tab. [5]) 
and radiative transfer (Tab. [7]) are in accord. The discrepancy in luminosities for 7 is due 
to lack of an outer envelope in the radiative transfer model, or presence of a Class source 
contaminating the 70 /im flux of 7, as discussed above. The discrepancy in luminosity for e 
is due to lack of an outer heated envelope in the radiative trans fer model. The accreti on of 
the envelope is not included as a source of heating (see §2.2.4 of iRobitaille et al.ll2006l ). 
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Fig. 8. — Spectrum of IC1 396Aa together with its best-fitting radiative transfer model from 
among those computed by iRobitaille et aLl (120061 ). The IRAC data are symbols with error 
bars; the IRS spectrum is the solid, jagged lines. Model curves are scaled downward by 10% 
for clarity. The model sum is a solid, smooth line. Individual model components are the 
outer envelope (dashed line) that contributes in the far- infrared, the inner envelope or 'disk' 
(dotted line) that contributes in the mid- infrared flux, scattered light (dash-dot line), and 
direct light (dash-dot-dot-dot). 
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The accretion rates, which are probably the most interesting quantities as they tell 
us how rapidly the stars are growing, are in the same order of magnitude {10~^Mq yr~^). 
In both the empirical and radiative transfer models, the mid-infrared flux is dominated 
by the envelope, whose luminosity is determined by the accretion rate; thus the accretion 
rates are fairly well constrained by the observations. Nonetheless, there are some significant 
differences. In particular, the empirical model yields an accretion rate 9 times higher than 
the best-fitting radiative transfer model for e. The massive outer envelope contributes 80% 
of the emergent luminosity in the empirical model for source e. The luminosities are also in 
general agreement, as expected since the models were chosen to match the observations at 
wavelengths that span the brightest part of the spectral energy distribution. For source 7, 
the empirical model yields a significantly higher luminosity, due to the cold, outer component 
that was required to match the 70 fim flux. The Whitney models cannot match the 70 //m 
flux, hence the poor flt quality for 7. 



-27- 



4. Class II Sources 

4.1. Disk Properties from the mid-infrared spectra 

Figure M shows the spectra of four Class 11 sources. Each spectrum was dereddened 
using the line of sight extinc tion derived in §2.51 and a nominal Ry = 3.1 extinction curve 



( iWeingartner &: Draind l200ll ). The extinction correction is small in the mid-infrared: for 
Av = 5.4, the optical depth r = 0.14 at the shortest IRS wavelengths (5.5 /im), 0.38 at 
the peak of the interstellar silicate feature (9.5 fim) and 0.06 at long wavelength (30 /xm). 
But the line of sight extinction is high enough to block much of the optical and nearly all 
ultraviolet light from these young stars, making it impossible to assign spectral types using 
the stars' colors. 

The mid-infrared spectra in Figure [H] are characterized by a modest near-infrared bright- 
ness, a very flat continuum from 3-8 /xm, a bright silicate emission feature at 9-11 /im, and 
a continuum beyond 15 /xm with a range of brightnesses. Table [5] lists the colors and silicate 
feature strengths for the sources. The silicate feature ranges from very strong (twice the 
continuum) to weak (20% of continuum) to undetectable. 

The radial and vertical structure of the disks can be complicated during the time of 
planet-building due to ejection of material from the planet growth zones as well as collection 
of material near stable resonances. Of particular importance for understanding formation 
of the Earth is the presence or lack of solid material in the ~ 0.5-2 AU 'habitable' region. 
For all the Class 11 sources with signiflcant emission at 24 /im, the color temperature of the 
15-40 /im emission is around 200 K. There is no evidence of the presence of a colder disk. 
Given luminosities of L*, the color temperatures suggest that the emission arises within 1.9 

1 /2 

LJ AU of the star, corresponding to the zone of terrestrial planet formation. Emission also 
arises from 5-10 /im, requiring dust even closer to the star. 
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Fig. 9. — Spitzer/IRS spectra of Class II sources in IC 1396A. Each spectrum is in uLi, 

units, divided by its value at 8 fim and offset vertically by 0, 2.5, 5, 7.5, 10 Lq, for sources 
IC1396A:^, p, Tr37 11-2146, Tr37 11-2037, and 21364964+5722270, respectively. The solid 
line is a simple blackbody through the near-infrared photometry to illustrate the approximate 
flux of the photosphere. 
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Mid-infra red spectra o f large samples of T Tauri a nd Herbig AeBe star s have been 
presented by iFurlan et al.l (l2006l . IRS: Taurus- Auriga) , iKessler-Silacci et al.l (120061 . IRS: 
Chameleon, Lupus, and other nearby star-fo r ming regions), iHonda et al.l (l2006l . Subaru: 
Taurus, TW Hydra), and Isiciha-Aguilar et al.l J2OO7I . IRS: Tr 37 and NGC 7160). Our spec- 
tra of Class II sources in IC 1396A, selected for 8 /xm brightness, are generally similar to T 
Tauri stars in other star-forming regions. Silicate emission features are present in emission 
for all 5 stars in Figure O There are two evident morphologies of silicate features. The top 
3 spectra show bright, distinct 9-11 and 18-22 fim features similar to other T Tauri stars. 
These features are produced by small amorphous silicate grains in an optically-thin portion 
of the disk, e.g. its upper surface or a flared or warped region illuminated by the star. 

The bottom two stars in Figure [HI IC 1396A:0 and (3, have much weaker and broader 
features at 9-11 fim and no evident feature at 20 /xm. Both sources are bright in the IRAC 
and MIPS images; in the near-infrared (Fig. [2]), both have have nearby confusing sources. 
Structured nebulosity could contaminate the spectra near the PAH and important ionic 
lines. 6 may have associated scattered light in the optical (but not in our Palomar or IRAC 
images). We applied local background subtractions when extracting the spectra, to minimize 
any contamination from the nebula or neighboring sources. The silicate features for these 
stars have low amplitude, but there are two peaks that are evident at the long and short- 
wavelength edges, making the feature appear 'horned.' The 'horned' appearance is not an 
effect of line-of-sight extinction. The interstellar absorption feature peaks at 9.5 fim and is 
broad. If the extinction to IC 1396A:^^ were Ay > 10, then any silicate emission feature from 
the star would have its 9-10 fim portion depressed, generating a very "red" emission feature 
containing only the A > 10 portion. For IC1396A:0, the two 'horns' are at 9.1 and 11.3 
fim. The 11.3 fim horn is at the appropriate location for olivine, Mga;.Fei_a:Si04. We now 
consider two possible explanations for the shape of the silicate feature. 

First, we consider that the peaks in the feature are due to small grains of specific miner- 
als. The 9.1 yum horn is at the approximate location for a phyllosilicate like montmorillonite 
[(Na,Ca)o.33(Al,Mg)2(Si40io)(OH)2nH20] or for amorphous quartz [Si02]. We calculated 
absorption efficiencies for forsterite, montmorillonite, and quartz usin g a continuous dis- 



tribution of ellipsoids (CPE) and laboratory optical constants (from iFabian et al.l 12001 



Glotch. Rossman. fc AharansonI 120071 : iHenning &: Mutschkd 119971 . respectively). A reason- 
able fit to the spectrum includes small grains (< 2 /xm) of both crystalline forsterite and either 
montmorillonite or quartz, with comparable abundance, together with continuum emission 
at temperatures ranging from 700 to 150 K. Figure [10] compares the fit and data. The silicate 
feature can be explained exclusively with crystalline silicates, though the mineralogy of the 
larger grains cannot be distinguished. The crystalline features are emitted predominantly 
from the hottest part of the disk. If the luminosity is 40 Lq (as inferred from the radiative 
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transfer model fit discussed below), then the crystalline silicates are located about 1.6 AU 
from the star. The terrestrial-planet-building materials around the stars IC1396A:6' and jS 
could therefore be rather different from most low-mass stars, even in the same globule (cf. 
Tr37 11-2037 and the other stars in the top of Fig. [H] that are dominated by amorphous 
silicates) . 

Second, we con sider whether grain growth can explain the shape of the silicate feature. 
Honda et al.l (120061 ) showed that if the grains are glassy olivine, the weaker silicate features 
can be explained by ~ 2 fim grains and the stronger silicate features can be explained 
by ~ 0.2 /xm grains. To determine whether larger grains can explain the spectrum of IC 
■ 396A:^, we calculated absorption efficiencies for a range of grain sizes for glassy olivine 
(IDorschner et al.lll995l ). then optimized the disk and size to match the observed spectrum. 
Figure fTUT a) shows the result. As noted previously, the larger grains have a flatter and 
wider spectra feature, eventually becoming featureless (and joining the continuum). A peak 
in the absorption efficiency of glassy olivine occurs around 8.6 /im for particles with radii 
greater than 12 /zm. As seen in Figure [TOl the peak is at wavelengths too short to match the 
observed feature; however, the general shape is very similar to the observed one, and wider 
mineralogical search may yield a better match. 

In summary, the distinctive silicate feature (weak, wide, and 'horned') of IC 1396A:0 and 
(3 could be due to crystalline mineralogy (with quartz or phyllosilicates providing the 'blue' 
horn of the feature) or by grain growth (with larger grains providing th e width of the feature 
and a lso explaining its weakness) . For w ide samples of T Tauri stars, iKessler-Silacci et al. 



J2006I . their Fig. 9). lHonda et~al] J2006I . their Fig. 2), and iFurlan elaP J2006I . their Fig. 6 
versus their Figs. 3-5) showed that the disks with weaker silicate features also have 'flatter' 
silicate features, where flatness is deflned based on the ratio of the amplitude of the 11.3 /im 
olivine peak relative to the 9.8 /im amorphous silicate peak. Our results in Figure M also 
show this trend. 



Furlan et al.l (120061 ) propose an evolutionary sequence based on the mid-infrared spec- 



tral morphology, driven by dust settling and growth. In this scheme, IC 1396A:6' and (3 
would belong to their 'Group D', the group with the greatest amount of grain growth. 
Oth er stars with similar spect ral shapes in th e star cluster include Tr 37 11-2006 and 24- 
515 (ISicilia-Aguilar et al.ll2007l ). More recently, iKessler-Silacci et al.l (120071 ) showed that the 
'flatness' of the silicate feature is also correlated with stellar luminosity (i.e. flatter silicate 
features tend to be present around more luminous stars). Our data also support this trend. 
In particular IC1396A:^ and f3 are relatively luminous stars (approximate types A2 and FO) 
with flat features, and the later-type stars hke Tr 37 11-2146 and 21364964+5722270 (K6 
and G5) have stronger features. 
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Another possibility, as discussed above, is that the mineralogy is different, with in- 
creased abundance of phyllosilicates or quartz around the higher luminosity stars. Phyl- 
losilicates could arise from the catastrophic disruption of a young asteroid, since aqueously 
altered silicates are common compositions of asteroids in particular in the outer main belt 
( jRivkin et al.ll2002l ). Amorphous Si02 may be produced from cat astrophic collisio ns as well, 
perhaps among protoplanets as suggested recently for one star (IRhee et al.ll2008l ). It must 
be acknowledged that low-sensitivity infrared spectroscopy, such as presented in this paper, 
does not uniquely provide independent signatures of grain growth or mineralogy, because of 
the wide range of possible size distributions and compositions. High-quality spectroscopy of 
a wider sample of T Tauri stars will shed more light on this issue, which is important since 
the a disk bearing phyllosilicates versus a disk with moderate amounts of grain growth are 
such suggest distinctly different evolutionary state for the disk. 



4.2. Properties of the visible spectra 



The optical spectra of the Class II sample show bright Ha emission lines, making them 
'classical' T Tauri stars. Figure [TT] shows so me of the spectra. The stars appear similar to 
those observed by (ISicilia-Aguilar et al.ll2005l ). suggesting they are part of the same popula- 
tion of T Tauri stars widespread throughout the IC 1396 region, of which IC 1396A is but a 
small part. Therefore, we do not elaborate upon their properties here. 



4.3. Association with Whitney protostar models 



The best-matching radiative transfer model from lRobitaille et al.l (120061 ) and lWhitney et al 



( 120031 ) was identified using the same procedure as used for the Class I sources ( §3.3.2p . As 
was found for the Class I sources, the best-matching models for the Class II sources also 
had fairly high x^, so the models are only to be used as a guide to plausible properties of 
the stars. Most of the mid-infrared-selected Class II sources are of order solar mass, with 
moderate-mass and slowly growing disks and negligible envelope. The two known T Tauri 
stars, LkHa 349a and c are fit by 2.2 and 0.7 Mq stars. The main difference between them 
is that LkHa 349c has a massive disk. 

One Class II source whose model stands out from the others in terms of luminosity and 
mass is IC1396A:^^. Its mid-infrared spectrum is very flat, and it could not be fit by low- 
mass models as could the others. The best fit model was a pre-main-sequence A star (SMq). 
Such a star has a large predicted photospheric decrement at 4000 Athat is not observed 
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Fig. 10. — Spectrum of IC 1396A:^^ together with three mineralogical fits. Panel (a) shows 
a combination of small (radius < 1 fim) crystalline olivine (forsterite; labeled 'fors') and 
phyllosilicate (montmorillonite; labeled 'phyll') particles. Panel (b) shows a combination of 
small crystalline olivine and large (radius 15 /xm) glassy olivine particles (labeled 'large oF). 
The minerals were placed at the same temperature, optimized in the range 500-800 K (as 
labeled) and combined with blackbody continua at 1200 and 120 K, as well as the extincted 
stellar photosphere. Individual components of the spectral model are shown as dashed lines 
with labels, and the model sum is shown as a solid line. 
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Fig. 11. — Optical spectra of Class II sources in IC1396A. Each spectrum was normalized 
by its median, scaled by 0.8, and shifted by unity from the previous for clarity. The bright 
Ha lines are evident for all sources, except LkHo; 349a, which has a P Cygni-like line profile 
as had been seen by Cohen & Kuhi (1979). 
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in the optical spectrum, and the predicted V-J=-0.15 whereas the observed V-J=3. These 
could be reconciled if the extinction Ay > 5 to this star; the best-fitting radiative transfer 
model has Ay = 8. The determination of the stellar mass from the infrared luminosity is of 
course very mo del- dependent; a spectral type in the A-F range are possible. There are other 
intermediate mass stars in Tr 37, including B stars in the IC 1396A region with surrounding 
nebulosity, so it is not implausible that IC1396A:6' is a young, intermediate mass embedded 
within a dense clump of material that is not its envelo pe. These properties plus the bright 



Ha line make IC1396A:^ plausibly a Herbig AeBe star (jWaters fc Waelkenslll998l ). 



The spectral energy distribution of IC 1396A:/5 is also relatively flat in the mid-infrared, 
but its 24 nm flux is not high, which led us to assign a Class II photometric type. The 
spectra presented in this paper confirm that /5 is a T Tauri star and not a Class I protostar. 
The best-fitting radiative transfer model has a large line of sight extinction. Ay = 6.8, 
explaining the lack of bright optical counterpart and infrared-dominated SED. The spectrum 
of /3 s hows faint PAH fe atures, similar to those seen around other young intermediate-mass 



stars (jSloan et al.ll2005l ). though bright nebular emission in IC 1396A may contaminate the 



spectrum at the PAH wavelengths. 

The spectral observations of Class II sources show that the mid-infrared-selected sample, 
identified in Paper I, is dominated by highly reddenned, pre-main sequence stars of a wide 
range of masses from 0.4-5 Mq and luminosities from 1-40 Lq. 
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Interaction between the Protostars and Globule 



5.1. Physical conditions in the globule 



The present configuration of the globule comprises an ionized layer on the side facing 
the O star (evident in Ha emission), a dense bright rim (evident in PAH emission), and a 
high-opacity region (evident in extinction and in mm- wave molecular emission). The dense 
region comprises a round 'head' (centered on LkHa 349) with a hole in its center, and 'trunk' 
extending to the SW of the head. 

Let us denote pr operties of the io nized layer with subscript i. From radio observations, 
its emission measure (IMatthewsl Il979l ) 



EM,, 



2000 cm"^ pc. 



(6) 



where Li ~ 0.07 pc is the thickness of the ionized layer; hence we infer the density Ui ~ 160 
cm~^ and, assuming temperature T, = 10"^ K, pressure Pi/k 



2 X 10^ cm ^ K. The ionized 



layer is maintained by the O 6.5 star, with the rate of ionizing photons incident on the 
surface of the globule equal to the case B recombination rate for gas with density n,. 

Let us denote properties of the globule interior with subscript g. The radius of the head 



Rg = 0.51 pc. From CO and ^'^CO measureme nts, the column den sity 



1.5 X 10^2 cm-2, 



and the velocity dispersion 6vg = 2.1 km s ^ flNakano et al.lll989l ). from which we infer the 



density Ug ~ 5000 cm ^ and pressure Pg = 3 x 10^ cm ^ K. Thus at present, the bright rim 
of the globule is in pressure balance with the ionizing layer. 

However, when the O star first generated its H II region and illuminated the globule, the 
pressure would have been much higher than that of the ambient gas. The initial conditions of 
the medium surrounding the O star are unknown, but th ey need not have been exc eptional 
in order to lead to the present configuration. For example, iGritschneder et al.l (120061 ) showed 
that an initially randomly turbulent medium evolves into a 'pillar' type morphology, like IC 
1396A or the famous pillars of M 16, due to the pressure and ionization of the O star. The 
region that became IC 1396A would have had an initial density higher than average. 

After the O star was born (~ 4 Myr ago), the pressure of the H II region drove a 
shock into t he proto-globule. Using the density of the ionized layer and globule now, and 
eq. 12-15 of lSpitzer. L.I (Il978l ). the shock speed would be ~ 2 km s~^. Such a shock would 
reach the center of the globule in 3 x 10^ yr. Of course the density of the proto-globule was 
initially lower, and the globule was initially larger. These effects approximately cancel: the 
shock velocity oc ^^'^ , while for a globule of constant mass the radius oc Uq so the shock 

1 /6 

crossing time oc , a very weak dependence on the initial density. When the shock reached 
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the center of the globu le, a strong over-de n sity d eveloped there. Numerical simulations, e.g. 
the landmark ones by iLefioch &: Lazarefl (119941 ). show both the central concentration and 
instabilities leading to smaller fragments within the globule. During the ~ 10^ yr period of 
compression, much of the globule mass becomes ablated in a 'cometary' phase. IC 1396A is 
presently in this phase, with stars of a wide range of mass forming from the over-pressured 
globule material. 



5.2. Outflows and Globule Destruction from the Inside 

The central hole in the head of the IC 1396A globule may be only one of several outflow- 
related structures in the globule. In the following subsections we discuss evidence for current 
outflows in the globule and for the shaping of the globule by outflows. 



5.2.1. Central cavity and LkHa 349 



In optical images, a relatively transparent hole, through which background stars are 
visible, surrounds the young star centrally located in the head of the IC 1396A globule. 
Some wispy nebulosity is present, with the appearance of refle ction nebulosit y . In m olecular 
line images, the hole appears relatively vacant of material. iNakano et al.l (119891 ) showed 
using high-resolution CO data that the central hole is indeed nearly devoid of molecular gas. 
In the Spitzer images, the hole also appears to be devoid of material, and furthermore shows 
a bright rim, with a particularly well-define 'wall' on its W side. The hole size is 45" x 58" 
(EWxN S). The central star, LkHa 349a, has a spect ral type F9e and effective temperature 
6165 K (j Cohen fc Kuhilll979l : iHernandez et al.ll2004l ). The extinction toward the star is is 
nontrivial (despite the nearly empty appearance of the hole). The extinction is mostly local 
to IC 1396A, as opposed to intervening line-of-sight interstellar dust (for which one expects 
less than 1 mag extinction at the distance of IC 1396), so the extinction curve may have 
high Rv ~ 5 typical of stars behind molecular clou ds (as opposed to the ra ndom ISM value 
of 3.1). Using the Ry = 5.0 (3.1) extinction curve, IHernandez et al.l (120041 ) find LkHa 349a 
to be a 5.4 (3.0) Mq star with luminosity 270 (45) Lq. The preferred, higher-luminosity 
solution is near the birth-line (age ~ 0.3 Myr) and consistent with the other indications of 
the star's youth such as the Ha line, which has a P Cygni profile indicating a wind with 
Vw > 300 km s~^. The width and shape of the emission and absorption lines suggest a very 
high rotation rate vsini ~ 190 km s~^, near the breakup speed; this may be related to its 
apparently high-speed wind. We note that LkHa 349a will be a late B - to early A- type sta r 
on the main sequence and is likely a very young Herbig AeBe star (iHerbig &: Radll974j ). 
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many of which have similar properties (reflection nebulos ity, self-absorbed H a winds) and 
are associated with star clusters with age of order 1 Myr (IHabart et al.ll2003l ). 



The location of the LkHa 349 inside the hole, and its high luminosity and youth, lead 
one to suspect the star could have cleared the hole. This can readily be shown to be plausible. 
If the star produced, for a time t,„, a niechan ical luminosity, L^, leading to a bubble of radius 
Ry,, then fICastor. McCrav. fc Weaverlll97,^ 



R^ = O-ULy^'rig^Hlr^ pc. 



(7) 



where t^^s = t^/lO^ yr, = n^/lO^ cm~^, and is in solar luminosities. If the stellar 
wind generated the presently-observed hole, i?^, = Rh = 0.09 pc, then we can solve for the 
mechanical luminosity 



.5000 cm-3 

and, for a wind speed v^,, the required mass- loss rate 



1 X iQ-hl, 



300 km s^^ 



5000 cm- 



MqJI-\ 



(9) 



Since Vw > 300 km s~^, the mass-loss rate required for LkHa 349a to blow a bubble in the 
globule is plausible, by comparison to oth er young stars. Class I sources with luminosity 
~ 10^ Lq have winds with ~ 0.1-1 Lq (jStahler fc Pallal 120041 ). Such a wind would only 
need to be active for ~ 10^ yr to provide momentum to sweep the cavity clear. The SED 
models for LkHa 349a and c indicate ages of 0.4 and 1.2 Myr, respectively, al lowing time for 



their wind phase to have passed. LkHa 349a is rotating exceptionally rapidly (IHessman et al. 



19951 ). and has little or no disk, perhaps enabling a relatively spherical and powerful wind. 
The star appears close to the stellar birth-line but is clearly exposed, suggesting that while 
young, the star is likely to have completed its main 'wind-bearing' Class and I phases and 
resides within the bubble that it blew. 

The existence of a bubble blown LkHa 349a provides supporting, circumstantial evidence 
that it could have been triggered by radiative driven implosion. If LkHa 349a and its bubble 
had existed before the formation of the O 6.5 star, the bubble would not have survived to this 
day. The shock front driven by the pressure of the expanding H II region would have reached 
the center of the globule long ago, as discussed above. Such a large-scale shock would likely 
have significantly distorted the cavity sent it downwind of the young LkHa 349a. 
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5.2.2. Water maser near 101396:^ 

Figure [1] shows water maser emission toward source IC 1396A:7. 22 GHz H2O masers 
are associated with Cl ass protostars, an d only rarely (4%) with Class I sources, in a 



survey of such sources (iFuruya et al.ll200ll ). The mid- infrared spectrum of IC1396A:7 has 
the characteristics of a Class I protostar with a rapidly-forming star accreting material from 
within a dense envelope. Is the H2O maser from an outflow caught 'in transition' from Class 
to I? Or is the maser spot from a clump of gas far from the protostar? 



As noted by IValdettaro et al.l (120081 ). the maser spot is significantly offset from the IC 
1396A:7 by 5000 ± 200 AU. The maser peak coincides with a faint infrared source in the 
IRAC images. The coordinates of the IRAC source are 2l''36'"07.8" 57°26'42.8", within 1" 
of the VLA maser position; we will refer to this source for convenience as IC 1396A:75. 
The fluxes of 7,, at 3.6, 4.5, 5.8, and 8 /xm are < 0.2, 0.6 ± 0.1, 1.5 ± 0.2, and 1.2 ± 0.3 
mJy, respectively; with most of the uncertainty due to blending with the much brighter 
source 7. The corresponding magnitudes are [3.6]>15.3, [4.5]=13.7, [5.8]=12.2, [8.0]=11.8. 
No counterpart was evident in the near-infrared images. By chance, the orientation of one 
of the mid-infrared (IRS) spectrograph slits covered the faint source. The spectrum of 7^ 
cannot be cleanly separated from that of 7, but it is clear that the emission from the source 
is not line-dominated. If the infrared source were a Herbig-Haro type object, where the jet is 
impacting a dense clump, then we would expect strong emission lines of II2. The IRS short- 
low spectrum (where the faint source can be separated from 7) covers 5.2-14 /im and would 
contain the pure rotational H2 S(2) through S(7) lines if they were present. Furthermore, 
there is no corresponding spot at the H2O maser position in the near-infrared rovibrational 
H2V=l-0 S(l) image. Therefore, the IRAC source corresponding to the maser is not a shocked 
clump of gas. It is also unlikely to be a background galaxy, based on the lack of PAH features 
that would be prominent at 6.2 and 7.7 /xm modulo some small redshift. The IRS spectrum 
of the faint source is dominated by continuum emission, and it looks roughly similar to that 
of the 7; the difficulty of separating the sources makes it impossible to study the spectrum 
in detail from the available data. 

Therefore, it appears that the H2O maser is not directly associated with the luminous 
infrared source IC 1396A:7, nor is it an interaction spot between the outflow and nearby 
molecular gas. The faint mid-infrared source IC 1396A:7;, may itself be a protostar powering 
an outflow responsible for the H2O maser. At 24 yum, a very bright source clearly corresponds 
in position to the luminous source 7. There is an extension from this source at 24 /im that 
corresponds to the faint source seen with IRAC, but it is severely blended: the offset is 7" 
while the beam size is 5". We scaled and subtracted the point spread function to estimate the 
24 /im flux of the fainter source as ~ 10 mJy, corresponding to [24]~7.1; a very conservative 
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limit is [24] > 5.9. Combining the mid-infrared colors of 7^ and comparing to Figure 3 of 
Paper I (i.e. comparing to the locations of Taurus and IC 1396A young stellar objects), 
the H2O maser source 7f, appears to be a low-luminosity Class I protostar. This brings 
back the original question of why a Class I protostar has an associated maser. Perhaps this 
low-luminosity source is a Class I/O transition object. 

The remarkable result is that the brightest (only known) H2O maser in the globule is 
close to the second-most-luminous {Lmir = I.6L0) Class I protostar in the globule, but it 
is actually positionally coincident with a separate, low-luminosity (O.O2L0) infrared source. 
The separation of 5000 AU is somewhat wider than a binary pair; if 7 has a mass (star 
plus envelope) of ~ IM© then the orbital speed would be 0.5 km s~^, comparable to the 
velocity dispersion of the globule. Class I sources have been found to have a wide distribu- 



tion of separations, and this potential pair would be among the widest (cf. IConnelley et al. 



20081 ) . There are other low-luminosity infrared sources in this field, so the total popula- 



tion of protostars may be much larger than the bright mid-infrared sources identified by 
eye in Paper I. Such low-luminosity objects have been found in other nearby star-forming 



regio ns that were surveyed by Spitzer , including tho se in dense cores in L1014 (jYoung et al. 



2QQ4 \. L1448 (IKauffmann et all 120051 ). and L1521F (IBourke et all 120061 ): they have been in- 
terpreted as sub-stellar protostars though the envelopes may contain more than IMq. The 
proto type of the class, L 1014-IRS, was searched for H2O masers on 7 epochs with no detec- 
tion ( jShirley et al.l 120071 ). The detection of an H2O maser from IC 1396A76 is the first for a 
low-luminosity object of this class. Masers luminosities hav e been found to correlate roughly 



with the bolo metric luminosity of their driving protostars ( iFuruya et al.ll200ll : IShirley et al. 



20071 ). with LH2O = 3 X 10 ^Lq and over an order of magnitude of scatter. One could 



simply use this scaling to estimate or constrain the luminosity of the source, despite the 
large scatter and wide range of variability (a factor of 180). For IC 1396A:7ft, we would 
obtain Lt,oi ~ O.l-IL© depending which epochs are averaged together. Such a luminosity is 
not inconsistent with the observed mid-infrared luminosity, though it does suggest that the 
much of the luminosity emerges in the far-infrared. 



5.2.3. Molecular Outflow from IC 1396A.-f 



The Palomar 2.12 fim H2I-O S(l) continuum-subtracted image shows two knots of emis- 
sion near IC 1396A:7. Figure [12] shows that these knots are located diametrically oppo- 
site each other along a line that passes through 7, a configuration strongly suggestive of 
a jet, powered by 7, that is striking molecular knots forming dense Herbig-Haro objects 
(IReipurth and Ballyl 120071 ). The short length of the outflow (in particular compared to pc- 
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scale outflows) is likely due to its propagation into dense walls of the globule. That the 
outflow is seen in H2 but not other tracers is also likely due to the dense material the jet 
interacts with. 

Further evidence for a molecular outflow is found from observations this source using 
the Arizona Radio Astronomy Observatory Submillimeter Telescope in the CO (3-2) line (23" 
beam, including both 7 and 7^. Figure [13] shows a prominent blue wing; a two-component 
flt of narrow line typical of the globule plus a second Gaussian with a full width at half 
maximum of 7.9 km s~^ and peak antenna temperature 3.7 K. 
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The presence of outflows from young ste ll ar ob jects in IC 1396 is not unique to the 
globule studied in this paper. iReipurth et al.l (120031 ) found HH objects and a 0.6-pc scale 
outflow in IC 1396N. The find that HH 777 emerges from the dense globule into the H II 
region. In the following section we discuss the potential influence of outflows on IC 1396A. 



5.3. Compartments within the globule due to outflows 

The internal structure of the globule around IC 1396:7 (Figure [T2|) is complicated, but 
it can be described as follows. There is a bright 'wall' located directly ~ 10" upstream (i.e. 
toward the central O star of the H II region) of 7, with a curved shape that appears to merge 
with two east-west fllaments that pass ~ 30" north and south of 7; the fllaments continue 
downstream of 7 for ~ 180' before merging with other structures in the globule. 

Using the equations above, assuming the upstream wall is a boundary between the wind 
and the turbulent pressure of the globule (or the ionized layer surrounding the globule, which 
has comparable pressure), the required mechanical luminosity of the wind ~ O.O2L0 and the 
mass-loss rate ~ 1O~^M0 yr~^ (assuming an outflow speed ~ 30 km s~^ typical of such 
outflows) if the wind has been active for 10^ yr. Such a luminosity is plausible, as 7 is a 
class I/O so urce with luminosity ~ 3OL0, for which outflow mechanical luminosities ~ O.ILq 



are typical (IStahler &: Fallal l2004j ). The shape of the cavity around 7, if it is wind-blown, 
suggests a signiflcantly lower density on its downstream side. If we naively use the wind 
bubble radius equation for the upstream and downstream directions separately, then with a 
size some 18 times longer in the downstream direction, the density would be lower by a factor 
of 10^, meaning the region downstream from 7 was already essentially empty, with n < 1 
cm~^, compared to the rest of the globule. Thus it seems unlikely that the structure around 
7 is an asymmetric, wind-blown bubble. Rather, its appearance and physical properties 
suggest a wind blowing into a high-pressure front of dense globule material upstream, and 
expanding relatively freely downstream into a region of much lower density. 

Could this shape be a wind bubble blown by 7 and other nearby protostars? We can use 
the same type of estimate as in the previous section to estimate the wind properties. The 
globule is very optically dark near 7 indicating a high column density. The H2O maser line 
is narrow, but we cannot use this as an estimate of the wind speed because masers are highly 
localized spots within the outflow. Instead, a lower limit to the molecular outflow speed is 
obtained from the range of velocities observed for the H2O maser, which spanned 5 km s~^ 
over the observed epochs (Fig. [1]). The actual outflow velocity is likely to be much higher 
depending on projection effects (as the masers are only individual spots in the outflow and 
can only under-represent the true range of velocities) and the bias that the H2O emission 
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Fig. 12. — Spitzer/IRAC color image of the region around the protostar IC 1396A:7. The 
3.6 /zm image is coded as blue, 4.5 /im as green, 5.8 /xm as orange, and 8 /xm as red. The 
Palomar/WIRC H2 2.12 /zm image is coded as light green. Protostars IC 1396A:7, 5 and e 
are labeled. Source 7^ is coincident with the 22 GHz H2O maser. A small 'x' indicates an 
IRAC image artifact exactly 4.88" (4 IRAC pixels) SE of 7. Dashed lines indicate the two 
shocked H2 blobs on either side of source 7. The scale bar is 30"=0.11 pc=0.36 l.y.=22,000 
AU long. 
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Fig. 13. — CO (3-2) spectrum of the protostar IC 1396A:7. 
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traces only dense, warm molecular gas and not the likely ionized and much faster material 
closer to the star. 

This empirical reasoning can be further extended to the rest of the globule. There are 
5 prominent 'voids' within the globule and its main SE extension. All of the bright Class 
I sources for which we obtained IRS spectra are located within voids. The one containing 
7 was just discussed in the previous paragraph. The one containing 6 and e is similar in 
many ways to the one containing 7, in size, shape, and orientation. Since 6 and e are so 
close, it is possible that their winds combine to form this void. Three large voids are located 
within the globule head. These include the one containing Lka that was discussed above, as 
well as one that contains protostar a. The third contains sources /i and (, with its E edge 
overlapping the globule edge; however, its NE boundary is clearly not the globule edge and 
has a curvature consistent with a roughly parabolic surface with fj, at its focus and C near 
its edge. 

Since the major voids all contain either mid-infrared-bright. Class I protostars, or the 
Herbig AeBe star in the globule hole, we suggest that the globule is being reshaped from 
the inside, with each protostar residing in a 'compartment' that it has blown via its outflow. 
Turning the argument around, we should see what influence each of the photometrically- 
defined Class I sources may have had on its envir onment. For this we consider the sources in 



the Class I region of the color-color diagram from lReach et al.l (120041 ). In order of 8 flux, 
these are a, 7, 6, e, r], (, /i. All except r] are accounted for in voids. Since r] is located in a 
relatively busy portion of the globule, it would be difficult to separate any potential cavity 
from other overlapping structures. There are other, poorly-defined or smaller voids within 
the globule. One such void in the SW 'trunk' extension is located just W of 6+e and just N 
of 7; this void is not as well defined and contains no mid-infrared source. Thus it appears 
that all (but one) of the Class I protostars are located within voids, and all (but one) of the 
voids contain Class I protostars or the Herbig AeBe star. 

The evidence suggests that protostellar outfiows are shaping the globule from the inside, 
while the pressure from the H II region maintains the bright rim in pressure equilibrium. 
The outfiows may have a double significance. Not only are they destroying the globule's 
structural integrity, leading to its further fragmentation and dissipation, but they may also 
drive and maintain turbulence, which wo uld otherwise dissipate quickly in the dense gas. 



within the globule (ILi fc Nakamurall2006l ). This induced turbulence may be what prevents 



the entire globule from collapsing into stars, limiting the star formation efficiency. 
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5.4. Proplyd-like objects near IC 1396A 

In addition to the protostars and wind signatures, the mid-infrared im ages contain 



featur es indicative of dense, starless cores. The term 'proplyd' was coined by lO'dell et al. 



(I1993I ) to refer to protoplanetary disks made observable as silhouettes in HST images of 
the Orion H II region: a generic description would include a central star, an elliptical, dark 
patch, and a bright rim. Without a central star, we have a proplyd-like object that could be 
starless or harbor a protostar. There are numerous dark globules wi thin IC 1396, ranging 



from the relatively large and bright ones like IC 1396A or IC 1396N (IReipurth et al.ll2003l ) 
to smaller fragments dispersed throughout the H II region. We mention here some of the 
smallest ones, because the large, bright-rimmed globules contain multiple embedded infrared 
sources and do not meet the spirit of the definition as being produced by the circumestellar 
material of a single young star. 

In the 10' region around IC 1396A, there are several small clouds that have distinct dark 
cores surrounded by bright rims and ofte n with cometary tails. These look similar to the 



proplyd-like objects in other H II regions. iDe Marco et al.l (120061 ) studied those in 8 regions 



finding that only one objects in M 17 co ntains a cent r al sta r and may be a real proplyd 



(joining those originally found in Orion). ISmith et al.l (120051 ) explains the paucity of such 



objects outside Orion as being due to the short lifetime of the proplyd phenomenon. (Only 
Orion contains numerous stars < 0.5 Myr very close to the O stars.) Are the proplyd-like 
objects in IC 1396 and other H II regions small, individual-star-forming globules, or are they 
transient molecular cloud fragments en route to photoevaporation? 

Figure [131 shows images of 5 proplyd-like objects in IC 1396, and Table [TUl tabulates their 
properties. For objects A, B, C, and E, the position angle of the vector from the object to 
the central O star (O* PA) is nearly equal to the position angle of the vector from the tail to 
the head (Tail PA), indicating the tails are material ablated from the head, as is the case on 
larger scale for the IC 1396A globule. The central holes are marginally resolved with IRAC 
at 8 /xm, having diameters 5.0 to 15" compared to the PSF diameter of 2.4". For each object 
there is a distinct, dark 'head' that could be due either to a lack of material (i.e. an empty 
hole) or extinction. If we attribute the 'heads' as being dark due to extinction, then we 
require an optical depth at 8 fim rg > 1, which corresponds to Ay > 10 and Nh ~ 1.8 x 10^^ 
cm~^. The heads are not resolved in the MIPS images; the 24 /im image is consistent in 
morphology with a suitably smoothed and resampled 8 /xm image (so the 24 /im morphology 
is consistent with either extinction or a central hole). The central volume density required in 
the heads are of order 10^ cm~^ and masses of order 0.1 Mq. The free-fall time is short (10^ 
yr) in such high-density gas, but the cores can be supported if they have turbulent velocities 
around 0.2 km s~^ or a significant magnetic field. The properties of the heads are typical 
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of small dense cores in nearby star-forming regions (IMvers. Linke. fc BensonI Il983l ). The 



analysis just presented is similar in nature to that of iDe Marco et al.l (120061 ) . who set the 
optical depth to Ha emission to of order unity and concluded their proplyd-like objects are 
low-density and transient. In contrast, we find significant optical depth at 8 /im, requiring a 
significantly higher density. The near balance between likely turbulent and magnetic pressure 
versus gravitational force suggests the cores may indeed be forming stars. Their dynamical 
configuration, with a central density enhancement, bright rim, and ablated tail, suggest that 
the central core density could be higher than indicated at the angular resolution presented 
here, in which case that portion of the dense core will almost certainly collapse into a small 
star. 

An alternative explanation for the proplyd-like objects is a standoff between a wind 
and the ionization front. In this case the deficit in brightness at the 'head' is due at least 
partially to a lack of material, instead of extinction. Using the same relations as in §5.2. H 
the required wind mechanical luminosity = 0.002-0.09 t~l Lq. An outflow power of this 
order is typical for protostellar sources with luminosity or order 1 Lq. If the luminosity 
of such a source were emitted predominantly at 24 (70) /im, its flux would be of ~ 0.5 
(1) Jy, which is clearly ruled out with the Spitzer data. The observed properties of the 
proplyd-like objects is therefore only consistent with being wind cavities if the embedded 
sources are significantly under- luminous (in photons) or high-powered (in wind). With the 
available information, we suggest the hypothesis previously mentioned, of an evaporating 
dense globule, is more likely. 

While we described the objects as 'proplyd-like' based on analogy between their mor- 
phology and those of the silhouette protoplanetary disks in the HST image of Orion, another 
analogy is the evaporating gaseous globules ('EGGs') that are seen in many star-forming 
regions (IHester et al.lll996l ). The primary physical difference between EGG and proplyd is 
the lack of a bright central star, although even some EGGs are found to have low-luminosity 
stars or even outflows when observed at very high resolution. A point of distinction between 
these EGG-like objects is that they are not connected to a larger 'pillar' or globule that is 
being eroded by the ionization front. Instead, the objects we present here are isolated. They 
were most likely small regions of relatively high density that are currently being illuminated 
by the star, just like miniature versions of the IC 1396A globule. The proplyd-like objects 
discussed here are sirn ilar to the Thackeray's globules in IC 2944, which appear to be starless 
(IReipurth et al.l 119971 ). That these small globules lack indications of star formation in the 
Spitzer images shows that the mechanism of triggered star formation may not operate on the 
smallest globules, which may be evaporating so quickly (perhaps due to their larger surface 
area to mass ratio) that stars to not have time to form in their cores. 
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Fig. 14. — Spitzer 8 /xm images of proplyd-like objects in IC 1396. Coordinates are give in 
Table [T0| panel (a) contains objects A and B; panel (b) contains object C; panel (c) contains 
object D; and panel (d) contains object E. An arrow in each panel shows a 30" long vector 
pointed toward the central O star HD 206267. All images are oriented with N up and E to 
the left. 
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6. Star formation in IC 1396A 
6.1. Pre-Main Sequence ages for infrared-bright stars near IC 1396 A 

To determine the a p;es of the pre-main sequence stars near IC 1396A, we used the models 



from ISiess et al.l (l2000l ) for a relatively high metallicity Z = 0.04, based on the assumption 
that the young Cep 0B2 region should be more enriched than the cloud from which the Sun 
formed some 5 Gyr ago. Figure [15] shows the location on a color-magn itude diagram of the 



l ow-m ass pre-main sequence stars within a 9' radius of IC 1396A from ISicilia-Aguilar et al. 



( I2OO5I ) and from this paper. The photometry and extinction correction are not accurately 
known, so the large spread in H-R diagram and the precise locations of individual stars 
should not be over-interpreted. Taken at face value, LkHac is very young (~ 0.5 Myr) and 
possibly coeval with the more-massive LkHa 349a (discussed in §5.2. ip . The location of 
the source IC 1396A:6' can be only roughly estimated since it is not detected at 5500A; for 
illustration we placed it at the location appropriate for V=R+.5=17.5. The mid-infrared 
sources tend to be younger than the optically-selected stars. Most of the stars are in the 1-5 
Myr range, including the infrared-selected Class II sources for which Palomar and Spitzer 
spectra were obtained. 



6.2. Generations of Star Formation in IC 1396A 



Apart from LkHa 349, which formed near the center of the globule, the triggered 
star fo rmation in smaller condensations may be expected to follow a sequence. For ex- 
ample, ISugitani. Tamura. fc Ogural (jl99ll ) found an age gradient such that T Tauri stars lie 
closer to the trigge r ing O star than IRAS sources, which tend to lie behind the bright rim. 
Matsuyanagi et al.l (120061 ) showed this in detail for an individual bright rim. Does a similar 
sequence obtain among the young, mid-infrared sources in IC 1396A? There is a handful 
of Class I sources within IC 1396A — more than the typically 1 IRAS source for bright rims 
where the sequential trend was first noted (ISugitani. Tamura. fc Ogural 1 19911 ). The youngest 
protostars in IC 1396A are likely those with the most 'extreme' infrared colors; these are e, 
6, 7, a, in order of increasing [8]-[24] (and [3.6]-[5.8]) color. The sources e, 6, and 7 are lo- 
cated just behind (0.02 pc) bright rims. On a larger scale, all of these protostars are located 
downstream of the center of the globule, and all but a are in the 'trunk' portion, known 
as SFO 36, which may have been swept back from the 'head' of the globule. The luminous 
Class I sources are further behind the bright rim than LkHa 349 and the hole it has carved 
from the center of the head of the globule. Thus it appears that the centroid location of the 
youngest protostars is significantly behind the main concentration of globule mass. 
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Fig. 15. — - Her tzprung-Russell diagram for IC 1396A, using theoretical models from 
(ISiess et al.l I2OOOI ) with Z = 0.04. Solid red lines are pre- main sequence Hayashi tracks 
for stars with masses (right to left) 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.8, 1, 1.3, 1.5, 1.9, 2.7, 5, 
and 7 Mq. The first 7 dashed black lines show isochrones for stars with ages (top to bot- 
tom) 0.01, 0.03, 0.1, 0.3, 1, 2, and 5 Myr; the lowest dashed line shows the zero-age main 
sequence (ZAMS). Two filled cyan circles linked by a cyan line labeled 'a' show the location 
of LkHa 349a for extinction curves with Ry = 3.1 and 5; the larger circle is the preferred, 
Rv = 5 result. A separate cyan filled circle labeled 'c' marks the location of LkHa 349c, 
on the 0.5 Myr isochrone and 0.6-0.8 Mq Hayashi track. Blue t riangles mark the location s 
of low-mass pre-main sequence stars within the IC 1396A field (jSicilia-Aguilar et al.l 120051 ). 
DodgerBlue-colored squares show the locations of the mid-infrared-selected (and optically 
detected) stars from Table (H IC 1396A:0 is labeled and includes a line for a 1 mag range of 
extinctions. Violet asterisks mark the locations of stars just upstream of the bright rim of 
IC 1396 A from Table [H 
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Inspired by this observation, we searched for Class II young stellar objects just upstream 
from the bright rim. Indeed, there is such a set of infrared sources. We measured the mid- 
infrared fluxes of several of the 24 /im sources, and f ound the corresponding e ntries in 
2MASS (J, H , K) and where possible in the tables by (ISicilia-Aguilar et al.ll2005l ) (optical 
magnitudes, spectral type) or the USN0-A2 catalog]^ Table [TT] lists the 4 sources with [8]jll. 
All of these source s are visibly exposed, despi te being mid-infrared selected. For the 3 stars 
with spectra from ISicilia-Aguilar et al.l (120051 ) , emission lines make them classical T Tauri 
stars, with ages > 10^ yr. A detailed study of such sources these sources is not possible 
without spectroscopy, because of variable extinction. One of the sources has essentially 
stellar color in the mid-infrared, but the optical counterpart is very faint (though it does 
have measured proper motion in the USNO catalog, indicating it is clearly not a background 
galaxy). 



Figure [T6] shows the distribution of sources photometrically identified as Class I and 
Class II using the Spitzer/IRAC and MIPS data from Paper I. The Class II sources are 
distributed throughout the region, while the Class I sources are clearly located within the 
globule. There may be a higher density of Class II sources toward the globule head and 
just upstream, though the number of sources is too small to determine if the over-density is 
significant. 

Note that we do not see the same trend of locations of near-infrared sources relative 



to mid-infrared that ISugitani. Tamura. fc Ogural (Il99ll ) found for some other bright rims. 
Whereas they found the sources with J — K > 1.2 are located preferentially upstream of 
the bright rim, we find instead that sources with J — K > 1.2, and also even redder stars 
with J — K > 1.8, tend to be located preferentially toward the globule. However, we do not 
believe this is a true association of large number of Class II YSOs with the densest part of 
the globule. Instead the distribution appears more consistent with a pattern of extinction 
of background sources, plus the Class I sources that are associated with the globule. To 
show this, we mapped the stars with more normal colors {J — K < 1.2), finding a fairly 
uniform distribution over nearly the entire region except toward the globule, which appears 
as a dramatic deficit of blue sources. In contrast, the distribution of red sources traces the 
globule accurately, and the number of red sources is exactly what is needed to 'fill in' the 
missing sources toward the globule. In §2.5l we described the extinction map generated from 
near-infrared source colors; this extinction map is very closely related to the distribution 
of CO and diffuse infrared emission. Similar extinction maps were obtained for IC 1396 
globules by iFroebrich et al.l (120051 ). Therefore, it would be challenging to use a simple near 



''From the USNOFS Image and Catalog Archive operated by the United States Naval Observatory, 
Flagstaff Station ( |http: / / www.nofs.navy.mil/data/fchpix/l ) . 
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Fig. 16. — Spitzer 8 fim image of IC 1396A with sources photometrically identified as Class 
I enclosed in red diamonds, and those identified as Class II enclosed in blue circles. 
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infrared color criterion to select near-infrared sources as T Tauri candidates for IC 1396A. 



With the available data, we can summarize that there are some mid-infrared-bright 
sources just upstream from the IC1396A bright rim, and most of these sources may be clas- 
sical T Tauri stars. A small minority of the upstream stars may be highly embedded sources, 
but they are much fainter and less numerous than the Class I sources in the globule. There 
is a much wider distribution of Class II sources throughout the entire, large IC 1396 region. 
Thus we see some evidence for an age gradient versus density: the protostars are associated 
with the dense gas and the Class II objects are dispersed. Two possible mechanisms have 
been suggested for age gradients of this sort, and both of these mechanisms can work in 
concert. Stars that form earliest are gradually exposed as the pillar recedes due to progress 
of the ionization effect and 'rocket' acceleration in response to the ionized evaporative flow 
toward the exciting O star. While this scenario does appear to match the configuration 
for IC 1396A, in this scenar io the young stars are expected to appear near the 'tips' of 
pillars (IHester fc Deschll2005l ). whereas we observe the presently-for ming stars further down- 



stream. The other scenario is sma ll-scale sequential star formation (ILefloch fc Lazarefllll994 



Sugitani. Tamura. fc Oguralll99ll ). where the earliest stars to form in a globule actually con- 
tribute to the collapse and triggering of the next generation. 



7. Conclusions 

We studied a mid-infrared-selected sample of young stellar objects from Spitzer imaging 
(Paper I) using spectroscopic observations from radio to optical wavelengths. The photo- 
metric classifications in Paper I are confirmed by the spectroscopy. The Class I sources 
have deep silicate absorption features and the vast majority of their luminosity arises in the 
infrared. None of them show H2O masers (characteristic of Class 0) or optical (Class II) 
counterparts. The only H2O maser associated with our original sample turns out to be a 
separate object, IC 1396A:7b, with much lower luminosity ~ 0.02Lq than the protostar we 
identified originally, IC 1396A: 7. The mid-infrared-selected Class I sources can be fitted 
with 0.1-2 Mq stars with ages < 2 x 10^ yr. 

Compared to stars forming in the more quiescent Taurus molecular cloud, the physical 
conditions of the globule from which the stars are forming include a higher radiation field 
due to the 06 star exciting the H II region and nearby young B stars. There is also a higher 
pressure due to compression by the ionization front from the 06 star and winds from the 
many young and still-forming stars. The wind from the ~ 5Mq star LkHa 349a has cleared 
a hole in the center of the dense globule, and outflows from the Spitzer-discovered Class I 
objects appear to have swept material into 'compartments' whose walls are over-pressured. 
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The spectroscopy results show that the Class I objects are brighter in the far-infrared than 
the theoretical models predict, suggesting they have a warmer envelope. Ice features in 
the Class I object envelopes differ somewhat from those seen around such objects in other 
star-forming regions. While the same set of features is present, their amplitude is lower, per 
unit silicate dust absorption, and the wavelengths are slightly shifted. Both of these effects 
suggest the Class I envelopes are warmer than in Class I sources in other star-forming regions, 
perhaps due to the warmer gas temperature in the star-forming material in the globule IC 
1396A. The somewhat different physical conditions of the Class I envelopes in this globule 
will affect the chemistry at the time when planet formation is just beginning. Because the 
Sun is thought to have formed near an O star, and most nearby star-forming regions lack an 
O star, astronomical searches for the physical conditions of the early solar nebula would be 
best conducted in regions like IC 1396. 

The stars photometrically classified as Class II are confirmed as such based on the 
presence of optical spectra with bright Ha emission lines and mid-infrared spectra with 
silicate features in emission. Most of the disk spectra are similar to other commonly-seen 
T Tauri star spectra, dominated by amorphous silicates. But the star IC 1396A:^^ has a 
particularly structured silicate feature with two 'horns' at 9.1 and 11.3 fim. The 11.3 fim 
horn can be explained by crystalline olivine. The origin of the 9.1 /xm horn cannot be 
uniquely identified, with reasonably good fits provided by quartz (Si02), a phyllosilicate 
(montmorillonite), or > 12 /im radius amorphous olivine grains. Similar spectra are evident 
in other wide samples of T Tauri stars. Such spectra show that the planet-building material 
around at least some stars is dominated by crystalline silicates and is significantly different 
from interstellar silicates. 

The star formation history in IC 1396A cannot be described as a single "burst" of star 
formation: Class I and II sources in and near the globule span a range of ages and masses in 
a small region (2 pc across). It appears that the earliest star formation in IC 1396 occurred 
approximately 5 Myr ago, based on the ages of the Class II sources and the central 06 star. 
Some stars ahead of the present-day bright rim of IC 1396 A have become exposed by the 
recession of the globule and by their own outflows. Recently, < 1 Myr ago, a powerful wind 
from the young AeBe star LkHa 349a blew a bubble through the center of the globule. The 
star formation in IC 1396A we observe now as Class I sources occurred within the last 0.2 
Myr, based on the ages of the Class I sources. 

The most recent star formation occurs preferentially behind the center of the globule, 
just inside bright rims. The specific trigger for the most recent star formation may have 
included not only pressurization of the globule by the O star exciting IC 1396, but also the 
pressure of outfiows and jets from the earlier generation of stars and in particular those with 
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powerful winds like LkHa 349a. The young stars within IC 1396A appear to be reshaping the 
globule from the inside, clearing 'compartments' within the globule. Gradually, the globule 
is being fragmented due to the exponentially increasing number of stars (if indeed they are 
sequentially triggered). The mass contained in walls of the compartments will eventually 
become small enough, and their overlapping structure chaotic enough, that further star 
formation in the globule will end. The limiting factor in the star formation efficiency of this 
one globule will be the dynamic effect of energy input from the stars formed within it. 
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Table 1. Source list and observing log 



Name*^ 


RA DEC 


Class^ 


IRS= 


DBSP'i 


GBT^ 


IC1396A:a 


21 36 46.6 


+57 29 38 


I 


SL+LL+H 





Y 


IC1396A:/5 


21 36 55.2 


+57 30 30 


II 


SL+LL 





Y 


IC1396A:7 


21"36 08.0 


+57 26 37 


I/O 


SL+LL+H 





Y 


IC1396A:5 


21 36 19.3 


+57 28 38 


I/O 


SL+LL 





Y 


IC1396A:e 


21"36 18.4'^ 


+57 28 30 


I/O 


SL+LL 





Y 


IC1396A:r7 


01/7 771 OS 

21 36 57.8 


+57 30 55 


I 


SL+LL 





Y 


IC1396A:^ 


"1 /l T7) 

21^36 39. 2^^ 


+57 29 53 


II 


SL+LL+H 


Y 


Y 


IC1396A:t 


"1 /? 77) OS 

21 36 36.8 


1 r" '^0 1 /oo// 

+57 31 32 


II 






Y 


IC1396A:k 


1 r) 771 r' /IS 

21 36 56.4 


+57 31 51 


II 






Y 


IC1396A:A 


1 /j 77) r" /I OS 

21"36 54.8 


1 r" ^0 ^ r\f r\r\H 

+57 30 00 


I 


SL+LL 





Y 


IC1396A:^ 


01 /j0^77)/^O S 

21"37 02.9 


+57 30 48 


I 






Y 


IC1396A:z/ 


21 36 57.9 


+57 29 11 


I/II 






Y 


IC1396A:^ 


2l''35'"57.9' 


+57°29'11" 


I/II 






Y 


IC1396A:C 


2l'*37'"02.3" 


+57°31'15" 


I 






Y 


LKHao4yc 


2l'^36'"49.5" 


+57°31'22" 


TT 

ii 


CT 1 T T 1 TU 


\r 
1 




LkHa349a 


2l''36'"50.7" 


+57°31'10" 


III 


SL 


Y 




Tr37 11-2146 


2l''36'"57.6" 


+57°27'32" 


II 


SL+LL 


Y 




Tr37 11-2037 


2l'^37'"07.0^ 


+57°27'01" 


II 


SL+LL 


Y 




Tr37 72-1427 


2l'^35™16.3^ 


+57°28'22" 


II 


SL+LL 


Y 




21361664+5728404 


2l'^36™16.6^ 


+57°28'40" 


II 


SL 


Y 




21364964+5722270 


2l''36™49.7^ 


+57°22'27" 


II 


SL+LL 


Y 




Tr37 72-875 


2l'^35'"49.8^ 


+57°24'04" 


II 


SL+LL 


Y 




21364398+5729287 


2l'^36'"43.9' 


+57°29'28" 


II 


SL+LL 


Y 




21362507+5727502 


21^36^25. F 


+57°27'50" 


II 


SL+LL 


Y 





'^Designations frorn IC13 96A:x from Paper I, Simbad: [RRY2004] x; designations from 
( ISicilia-Aguilar et al.ll2005l ). Simbad: [SHB2004] Trumpler 37x-y; numeric designations 
from 2MASS All-Sky Point Source Catalog 



'^Classification based on mid-infrared colors (Paper I). 
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^Observed with the Spitzer/Infrred Spectrograph; SL=5.2-14 /im, LL=14-35 /im, 
H=high-res (10-35 /im). 

'^Observed with the Palomar 5-m/Double Spectrograph; Y=good spectrum, 
o=observed but not detected (no optical counterpart). 

^Observed with the NRAO Green Bank Telescope for 22 GHz H2O masers. 
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Table 2. Brightness (magnitudes) of observed sources 



Name 


B 


V 


R 


J 




H 


K 


s 


[3.6] 


[4.5] 


[5.8] 


[8.0] 


[24] 


IV^ -LOiyKJiV.LX 








1 7 




15 


31 


± O 


51 


Q 858 


8 44 


7 1 Q 


fi 98 
















12 


47 


1 1 

± ± 






Q Q7 


8 75 


8 ns 


5 04 


TCI '^QfiA-^ 








1 9 


76 






1 3 


95 


10 56 


9 22 


7 84 


6 80 


2 24 


IC1396A:5 








> 


20 






15 


76 


12.66 


10.9 


9.51 


8.44 


3.73 


IC1396A:e 








> 


20 






16 


95 


12.57 


10.9 


9.16 


8.55 


4.27 


IC1396A:C 












16 


48 


13 


62 


10.84 


9.86 


9.02 


8.84 


4.48 


IC1396A:?7 








17 


08 


15 


41 


14 


05 


12.40 


11.36 


10.6 


9.15 


4.93 


IC1396A:A 








17 


89 


15 


72 


14 


06 


11.92 


11.55 


10.50 


9.98 


5.91 


IC1396A:^ 






17.0 


11 


73 


10 


21 


9 


50 


8.414 


7.85 


7.21 


6.53 


4.07 


LkHa349c 


17.0 




12.4 


10 


59 


10 


81 


10 


13 


9.600 


9.23 


8.77 


8.65 


4.94 


LkHa349a 


15.2 




11.4 


11 


71 


8 


95 


8 


62 


8.440 


8.44 


8.25 


8.14 




Tr37 11-2146 


17.1 


16.92 


15.6 


12 


07 


11 


52 


10 


82 


9.932 


9.69 


9.13 


8.34 


4.89 


Tr37 11-2037 


16.3 


16.03 


15.1 


12 


58 


11 


72 


11 


29 


10.87 


10.59 


10.2 


9.51 


4.49 


Tr37 72-1427 




18.08 




14 


03 


13 


07 


12 


67 


12.08 


11.89 


11.4 


10.8 


6.81 


21361664+5728404 








17 


02 


15 


06 


13 


98 


12.29 


11.57 


10.8 


10.0 


7.22 


21364964+5722270 


16.8 




14.9 


12 


84 


12 


07 


11 


79 


11.38 


11.51 


10.9 


10.3 


5.73 


Tr37 72-875 




18.37 




14 


17 


13 


04 


12 


51 


11.60 


11.31 


10.7 


10.0 


6.73 


21364398+5729287 








14 


34 


13 


07 


12 


38 


11.64 


11.24 


10.8 


10.6 


6.35 


21362507+5727502 








14 


84 


13 


94 


13 


39 


12.39 


11.96 


11.5 


10.4 


7.58 
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Table 3. Derived properties from optical spectra' 







Rd/Ra 






-^y 


A pre 


IC1396A:a 


1.12 


0.134 


3.6 


10.1 






IC1396A:;g 


0.86 


0.12 


4.1 


14.0 


FO^ 




IC1396A:7 


1.05 


0.19 


2.3 


(1.5) 






IC1396A:(5 


0.87 


0.16 


3.0 


4.2 






IC1396A:^ 


1.28 


0.071 


6.0 


3.8 


A2b 




LkHa 349c 


2.85 


0.092 


5.0 


3.6 


G9^ 




Tr37 11-2146 


4.71 


0.10 


3.4 


1.5 


K6.0 


0.9 


Tr37 11-2037 


11.9 


0.14 


3.3 


1.5 


K4.5 


2.5 


Tr37 72-1427 


4.8 


0.15 


3.1 


1.5 


Ml.O 


2.2 


21361664+5728404 


1.13 


0.177 


2.5 


3.6 






21364964+5722270 


5.02 


0.066 


6.2 


1.5 


05^ 




Tr37 72-875 


2.75 


0.145 


3.3 


2.4 


M0.5 


8.4 


21364398+5729287 


2.41 


0.071 


6.0 


7.1 






21362507+5727502 


3.26 


0.17 


2.7 


1.7 


MO.O 





^Ha and H/5 fluxes were measured for each source. For the Class I ob- 
jects, the spectra were extracted at the position of the infrared source and 
integrated over the appropriate wavelengths, but no optical counterparts 
were evident (so the lines are most likely nebular). Spectral types are from 



Siciha-Aguilar et all (j2005|) 



'^approximate spectral type based on luminosity 



Table 4. 22 GHz Water masers in IC 1396A 



Date 


ppeak 


Fh2Q 


Lh20 




(Jy) 


(Jy km s~^) 


(IO-^Lq) 


2006 Aug 07 


0.030 


0.02 


0.2 


2006 Sep 28 


0.042 


0.06 


0.7 


2006 Nov 18 


0.032 


0.02 


0.2 


2007 Jan 24 


0.29 


0.28 


3.7 
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Table 5. Mid- infrared spectral properties of protostars in IC 1396A^ 



Source [8]-[24] [3.6]-[5.8] as-s AFio/F AF15.2/F AFg.o/F 



Class I Sources 



IC1396A:q; 


4.27 


2.66 


-0.32 


-0.64 


-0.13 


-0.08 


IC1396A:7 


4.55 


2.72 


+0.19 


-0.70 


-0.10 


-0.06 


IC1396A:(5 


4.70 


3.15 


+0.16 


-0.79 


-0.18 




IC1396A:e 


4.27 


3.41 


-0.04 


-1.00 


-0.25 


-0.17 


IC1396A:77 


4.22 


1.75 


+0.24 


-0.31 






IC1396A:A 


3.57 


1.26 


-0.12 


-0.15 






Class II Sources 


IC1396A:/3 


2.99 


1.07 


-0.26 








IC1396A:^ 


2.45 


1.19 


-0.22 


0.208 






LkHQ;349a 


1.62 


0.82 


-0.30 








LkHQ;349c 




0.18 


-0.49 








Tr37 11-2146 


3.44 


0.79 


-0.38 


0.817 






Tr37 11-2037 


5.01 


0.63 


-0.31 


0.507 






Tr37 72-1427 


4.08 


0.58 


-0.31 








21364964+5722270 


4.61 


0.38 


-1.97 


1.94 






Tr37 72-875 


3.31 


0.87 


-0.03 








21364398+5729287 


4.32 


0.78 


-0.95 








21362507+5727502 


2.83 


0.88 


-0.15 


0.625 







^[a]-[b] is the color (magnitude difference) between bands a and b; cts-g is the spectral 
slope Fi, oc i/" from 5 to 8 /im; AF^/ F is the amplitude emission (or absorption, if 
negative) feature, centered on wavelength A (//m), relative to the local continuum. 
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Table 6. Ice Column Densities 



Source 


Af(H20) 


iV(C02) 


T9.7 




IQlS cm-2 


% H2O 




IC 1396A a 


1.6 (0.2) 


14 


1.12 (0.11) 


IC 1396A 7 


1.2 (0.2) 


34 


1.41 (0.14) 


IC 1396A 5 


4.3 (0.8) 


23 


1.96 (0.20) 


IC 1396A £ 


6.8 (1.7) 


26 


9.64 (1.00) 



Table 7. Empirical model for Class I protostars in IC 1396A 





IC1396A:7 


a 


6 


e 


Units 


Tsil 


1.7 


1.5 


2.1 


4.0 






24.5 


21.6 


30.2 


57.6 






0.5 


0.9 


2.0 


3.4 




^acc 


27 


24 


209 


279 


Rq 


R2 


3.3 


3.3 


6.8 


5.9 


AU 


Rs 


9.5 


1.8 


1.6 


1.7 


10=^ AU 


n2 


0.9 


0.8 


0.6 


1.2 


10^ cm~^ 


n3 


0.6 


6.4 


16.3 


24.8 


10^ cm-=^ 


M 


2.3 


2.0 


4.1 


7.3 


IO-^Mq yr-i 


Li 


2.7 


2.7 


0.6 


0.8 


Le 


Ltot 


12.3 


5.0 


1.8 


2.6 
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Table 8. Radiative transfer models for Class I protostars in IC 1396A 





IC1396A:7 


a 


5 


e 


V 


C 




Units 


ModeP 


3016360 


3000224 


3017673 


3016207 


3003269 


3013344 










i.4i 


Q 

O.D 


D.O 




i.o 






M 


91.4 


2.39 


2.55 


0.76 


19.3 


27.8 


10~ 


yr-i 


R* 


10.0 


6.13 


4.94 


4.10 


3.53 


3.39 


Rq 






1.88 


0.38 


0.15 


0.11 


0.20 


0.43 


Me 






4321 


3503 


2839 


2581 


3062 


3650 


K 




Age(*) 


11.1 


0.14 


0.16 


0.18 


4.00 


19.3 


10^ 


yr 


Ltot 


36.2 


5.37 


1.58 


0.67 


0.99 


2.68 


Le 




i 


69.5 


31.8 


41.4 


56.6 


41.4 


63.3 


o 




A int 


65.2 


62.9 


64.3 


73.6 


26.5 


30.4 


maj 


I (to central star) 


Av 


7.4 


1.4 


8.0 


7.4 


4.4 


4.4 


ma^ 


I (line of sight) 



^ Model identification from iRobitaille et al.l (120061 ). 



Table 9. Radiative transfer models for Class II protostars in IC 1396A 





21364964+5722270 


IC1396A:/3 


LkHa 349c 


LkHa 349a 


IC1396A:6I 


Tr37 11-2037 


Units 


Model'' 


3004629 


3015197 


3002727 


3005288 


3017958 


3007956 




xi 


3.9 


0.6 


8.7 


6.1 


5.6 


1.1 






0.009 


31.9 


25.9 


0.67 


0.00007 


0.26 


lO-^M© 




0.007 


68.2 


70.9 


0.01 


0.002 


1.53 




-R* 


1.77 


2.34 


4.41 


3.63 


2.67 


2.24 


Rq 


M* 


1.50 


0.80 


0.70 


2.28 


4.99 


0.80 


Mq 


T* 


5247 


4123 


3988 


4829 


16,450 


4124 


K 


Age(*) 


91.6 


10.1 


4.1 


12.4 


31.8 


11.7 


10^ yr 


Ltot 


2.1 


7.59 


4.73 


6.43 


46.9 


1.31 


i© 


i 


41.4 


69.5 


69.5 


81.4 


41.4 


18.2 


o 


J^int 


0.002 


0.016 


0.85 


0.003 


0.001 


0.011 


mag (to central star) 




2.8 


6.8 


1.80 


1.60 


8.0 


1.40 


mag (line of sight) 



''Model identification from IRobitaille et al.l 1I2OO6I) . 
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Tablc 10. Proplyd-likc objects in IC 1396 



Name 


RA 


Dec 


0* Dist 


0* PA 


Tail PA 


h 




Rhole 


Ltail 








(pc) 


n 


n 


(MJy sr 




(AU) 


(AU) 


A 


2l'^36"*25.2^ 


+57°39'48" 


5.0 


207 


196 




1.4 


2000 


10500 


B 


2l''36'"24.2" 


+57°39'22" 


5.0 


207 


199 




1.5 


1900 


9800 


C 


2l''38™13.7" 


+57°40'30" 


2.6 


242 


240 




2.4 


3300 


14000 


D 


2l''39™12.2" 


+57°24'07" 


1.2 


112 


35 




3.1 


5600 


74000 


E 




+57°26'59" 


4.7 


174 


161 




5.5 


1900 


25000 



Table 11. Selected Infrared Sources Upstream from IC 1396A 



2MASS [24] [8.0] [5.8] [4.5] [3.6] K, H J V 



21370936+5729483 6.70 9.83 10.31 10.84 11.19 11.84 12.33 13.39 18.2 

21371054+5731124 7.35 10.06 10.41 10.89 11.25 11.73 12.17 13.09 16.74 

21370649+5732316 6.99 9.91 10.44 10.69 11.19 11.91 12.33 13.32 17.85 

21371172+573329 6.95 10.78 10.73 10.88 11.06 10.90 11.12 11.94 17.4 

21370802+5734094 7.66 7.970 ... 8.01 7.74 7.86 8.37 9.62 14.5 



